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I. INTRODUCTION 


This report contains the results of a two and one-half year analytical 

* 

study of forced convection boiling of potassium. The study was undertaken in 
an effort to provide mathematical descriptions of the controlling processes in 
a forced-convection boiler tube, where such analyses were not already available. 
The analytical descriptions, when verified by comparison with experimental 
data, would then contribute to the rational engineering design of forced convection 
potassium boilers for space applications. 

All of the boilers presently under development or being proposed for use 
in Rankine-cycle engines for space auxiliary power are of the forced-convection 
once-through design. Such boilers receive subcooled or saturated liquid at the 
inlet of a continuous tube, and discharge dry, possibly superheated vapor from 
the outlet. In steady state operation, the flow past each point in the boiler is 
described by a fixed vapor quality. Typical flow regimes in such a boiler may 
be (in order, from inlet to outlet): 

1. All liquid. 

2. A bubbly, quasi-homogeneous mixture of liquid and vapor. 

3. Separated annular flow with liquid on the wall and vapor in the core. 
Liquid droplets may be entrained in the vapor. A wave structure may 
exist on the liquid film, and liquid fragments may be stripped from 
the waves by the vapor stream, ejected from the film by nucleate 
boiling, or deposited in the liquid layer from the vapor. 

4. Wetted liquid rivulets or droplets on the wall, wet or dry vapor in the 
core. 

5. Film boiling droplets . 

Under some conditions, film boiling droplets slide along the wall (on 
a film of their own vapor) while the vapor in which they are carried 
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superheats above the saturation temperature. Departures of the 
mixtures from thermodynamic equilibrium may be quite large. 

6. Fog flow. 

Under some conditions, fog flow, with a dispersion of fine droplets 
evaporating in superheated vapor, may occur. 

7. Vapor superheating. 

After all of the liquid is vaporized, dry vapor superheating, a process 
identical to the forced convection heating of any gas, occurs. 

This sequence of processes does not appear in every boiler. In some cases, 
flow regimes listed do not occur, or are replaced by others. A number of 
transition processes can occur which fall in between the flow regimes listed 
above. Some of these are: 

1. Surface boiling with the mixed mean of the liquid subcooled. 

This condition is readily obtained with low turbulence flows of poor 
thermal conductors such as organic liquids. It is probably not 
usually encountered in flows of liquid metals. 

2. Superheating of liquid above the saturation temperature without 
vapor generation. 

3. Transition boiling. 

At the end of the annular liquid film (or rivulet flow) , there is 
usually a sudden and very large drop in heat transfer conductance. 
However, some experiments have shown a relatively gradual decrease 
in conductance, indicating that liquid still partially covers the wall, 
possibly because spray is deposited from the vapor core, or wetted 
patches are swept downstream from the ruptured liquid layer. 

The literature of forced convection boiling includes many experimental 
investigations of n burnout" or "critical heat flux. " 

"Burnout” in such experiments does not occur uniformly over the length 
of the tube, but at a distinct axial location, usually at vapor qualities signifi- 
cantly below one hundred percent. It is clear that the term "burnout" has no 
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meaning in boilers of the type being discussed in this report since these boilers 
are required to produce dry vapor (one hundred percent quality) at the outlet 
without failure. The usual interpretation of the "burnout" phenomenon is that 
the annular film of liquid on the inner tube wall terminates at the vapor quality 
and axial location where the tube failure occurs. In the "once-through" forced 
convection boilers with which this report is concerned, this downstream termination 
of annular flow is an important location, and the parameters which determine it 
are very much of interest, but the boiler must operate around this phenomenon, 
and not be limited by it. 

In order to prevent entrainment of liquid into the vapor or to force the 
deposition of entrained liquid on the wall, boiler tube inserts are usually used. 

These may be of many configurations, but almost all rotate the fluid about the 
tube axis to promote phase separation by centrifugal force. The most commonly 
used inserts are: 

1. A plug with helical passages on the periphery. 

This has been used in the low quality region where the decreased 
flow area may be necessary to bring the vapor velocities up, 
thereby strengthening the separation forces. 

2. One or more helical vanes wrapped about a cylindrical centerbody. 

The centerbody lies along the tube axis, and may carry temperature 
and pressure instrumentation in experimental boilers, but its 
primary purpose is as a support for the vane. 

3 . A twisted tape helical divider . 

This divides the boiler tube into two passages. It requires no 

centerbody for support. 

4. A helical wire insert. 

This is simply a wire coil spring whose wire diameter and coil 
pitch can be chosen over wide limits . Sometimes a centerbody is 
used with wire coils to locate the coil axially (through struts) . If 
the centerbody is as large as the inside diameter of the coil, the 
action of the insert is in most 1 respects identical to that of a single 
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vane helix. 

5. Single vane helix without centerbody. 

Such an insert has been used but its fabrication and support in the 
boiler both present difficulties . 

It is clear that these inserts must have a significant effect upon the 
processes occurring in the boiler, and not all of the inserts affect a process in 
the same way. The nature of the insert, if any, must be considered in the 
analysis of forced convection boiling processes. 

Not all of the processes to be expected in a forced convection boiler are 
of equal importance to the designer. Prom a practical point of view, it 
probably makes little difference in the process of boiler design whether the flow 
early in the boiler is characterized by a homogeneous mixture of vapor and liquid, 
or by annular separated flow. In either case, the heat transfer conductances to 
be expected are higher than the already large liquid metal conductances of the 
inlet flow. (Which of these two regimes exists at a location may have a large 
effect on the amount of liquid entrained in the vapor at downstream locations, 
however.) On the other hand, it is of extreme importance to the designer to 
know what are the conditions controlling processes in the region of, and past the 
end of, the annular film, since the overall length of the boiler is strongly affected 
by the location of this region. 

In the study reported here, primary emphasis has been placed upon the 
analysis of processes which would be of maximum interest to a designer of 
practical equipment (although not all such processes have been considered) . 

As may be inferred from the foregoing discussion, most forced convection 
processes are, in detail, extremely complex, and a comprehensive analysis of a 
particular process might be expected to yield quite complex and unwieldy results 
and, in addition, be filled with assumed functional relationships of unknown merit. 
Such analyses have been avoided in the present work. The analyses presented here 


4 



are based on simple models of the flow, and the idealizations upon which the models 
are based are described in detail. The analytical results are simple enough in 
most cases so that the interrelationship between controlling parameters may be 
understood readily. Such simplicity, of course, limits the applicability of the 
analytical result to ranges in which the real flow is suitably approximated by 
the idealized model. In many cases, it is possible to establish the applicability 
range analytically, and this has been done. In other cases the reader must use 
his judgment as to the applicability of a given analysis for a particular set of 
parameter values. 


5 



II. ANALYTICAL STUDIES 


A. Flow Through Tubes Containing Helical Wire Inserts 

1) Helical Wires With Very Short Pitch-Tube Diameter Ratios 

The investigation of single phase pressure drop in tubes containing 
helical wire coils has shown that for very short pitch coils a good 
approximation of the pressure drop may be obtained by considering 
flow in a tube containing annular rings. 

A wire ring whose diameter is small compared to the diameter of 
the tube lies in the boundary layer of the flow in the tube. The 
boundary layer flow, which is separated from the tube wall down- 
stream of the wire may or may not reattach to the wall before the 
next wire is encountered, depending upon the helix pitch. For 
the longer pitches used in General Electric experiments described 
in References 1, 2, 3, and 4, the boundary layer reattaches to the 
wall before the next wire is encountered, but reattachment does 
not occur for the shorter pitches . 

In the case of the model of axial flow in a tube with annular rings, 
one would expect that a maximum would occur in the pressure drop 
versus ring spacing curve at the shortest pitch for which the flow 
still reattaches to the wall between rings, and that the pressure 
drop function would decrease somewhat with closer spacing. The 
pitch at which this maximum occurred would be dependent upon the 
relative diameter of the rings, since this height controls the 
reattachment length. 
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In order to relate the pressure drop in such a tube to other configurations 
which have been studied, a brief discussion will be made of the flow 
through a tube containing annular rings . 

Where the ridge spacing is wide compared to the height of the ridge, a 
relatively small stalled area in the tube wall boundary layer will exist 
upstream of the ridge, and a relatively large recirculating stall will lie 
downstream of the ridge. If the ridge spacing is wide enough, the main 
stream will expand behind the ridge and reattach to the wall before the 
next ridge is encountered. An extensive study of two-dimensional plane 
stalled regions is given by Abbott in Reference 5. Although the present 
case is axially symmetric, the flow patterns are similar. 

Downstream of the reattachment point a laminar boundary layer grows 
and changes to turbulent, just as in the entrance region of a tube, 
although the turbulence level is very high due to the eddies in the free 
shear layer coupling the main stream and the recirculation. It will 
be recognized that in the axially symmetric case the flow described is 
similar to that in the familiar orifice meter, and it will be remembered 
that the pressure fluctuations in the orifice meter are due to very 
high turbulence levels downstream of the orifice plate. 

Where the spacing of the ridges is short compared to their height, the 
flow will not reattach to the wall. Instead it will expand very little 
between ridges and the apparent tube diameter will be essentially the 
minimum diameter of the ridges. In the cavities between, a recirculating 
flow will be coupled to the main stream by a free shear layer, again 
of very high turbulence level. The pressure loss between ridges in 
this latter case will be less than in the case of reattachment, because 
there is no acceleration between the ridges. 


7 


Between the reattached flow case (wide spacing) and completely 
separated case (short spacing) , a region will be encountered in which 
the flow oscillates between these two cases with time, or in which 
some cavities are completely separated and others reattached, leading 
to pressure drop characteristics between the two extremes. 

Corrugated tubing and piping, which has been in use for many years 
and which has been investigated fairly extensively, has flow patterns 
quite similar to those just described, although it has only been fairly 
recently observed that the friction factors depend as strongly as they 
do upon corrugation pitch and height (reference 6) . The ranges of 
some of the early friction factor correlations for these tubes are very 
narrow. 

In the case of tubes containing short pitch helical wire inserts, the 
flow pattern is virtually the same, although now, since the recirculating 
regions are connected a pressure gradient will exist parallel to the 
wire, and the fluid in the stalled regions will drift helically downstream. 
Thus, in a flow visualization experiment using dyes, the concentration 
of dye in the stall downstream of the wire has a helical motion. Such 
a photograph, taken from Reference 3, is shown in Figure 1. The 
main flow need not have any significant rotation at all for this helical 
motion to be observed, since it is the local pressure gradient in the 
stall which is causing the motion. 

A comparison of a non-wetting two phase flow such as mercury-air 
in glass in a short pitch helical wire system and in a helical vane 
system of the same pitch and quality shows that the liquid component 
on the wall lies downstream of the helical wire, and upstream of the 
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Tube I.D. = 0.87S in. 

X = 0.84 

p/D = 0.83 

Wire Dia. = 0.125 in. 

A? = 0.348 psi/ft 

d w /D = 0.143 



Figure 1. Illustration of high quality flow on a wire coil insert. (Courtesy of 
General Electric Space and Missile Division. Cincinnati, Ohio.) 


helical vane. Secondary flows and the leakage between the vane and 
tube wall cause the liquid to pile up against the vane. This comparison 
also shows much less fluid against the wall in the helical wire system 
since the helical wire does not remove the liqu'id from the main stream 
as effectively. 


Abbott has given a relationship for the approximate length of the stall 
behind a single backward facing step, which may be used to compute 
the minimum pitch for which a reattached flow may be expected. 

Since Abbott's relationship is for plane rather than axially symmetric 
flow, and does not consider either the effect of the next downstream 
ridge or the contraction upstream of the ridge, the minimum pitch 
obtained from his work would be expected to be less than that actually 
encountered. The criterion for reattached flow based on Abbott's 



The line of points coinciding with this relationship is shown on the 
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friction data due to Sams and Converse in Figure 2. The data 
point represented by Figure 1, which clearly shows reattached 
flow, appears at p/D = .83. Also shown on this graph are several 
friction factor measurements taken from corrugated tube pressure 
drop experiments . ^ ^ 

The minimum pitch for a helical wire coil system is that in which the 
adjacent coils touch each other, i.e., p/D = d^/D. Two of the points 
due to Sams are for this case. It is clear that a large portion of the 
increase in friction factor for close-wrapped coils is due to the decreased 
flow area, since the effective diameter of the tube is smaller by twice 
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SMOOTH TUBE 






the thickness of the wires. It may be shown that the ratio of the 
friction factor based upon the tube diameter to that based upon the 
wire coil inside diameter is 



where 

£ , friction factor 

D , tube diameter 

D , coil inside diameter 
m 

d .wire diameter 
w 

The quantity , is a function also of d /D where d /D is 
D w w 

m 

essentially a measure of tube "roughness. " This quantity has been 
taken from Sams' data and corrugated tube data from References 6, 
8,9, and 10. The predicted curve of minimum pitch friction factor, 
based on these sources is also shown on Figure 2 as line I-I. 
Predicted asymptotes for some of the helical coil cases are indicated 
on this line. 

In the cases where the flow reattaches to the wall the pressure drop 
due to the annular rings is much greater than that due to wall friction 
and it may be expected that the pressure drop would b e directly 
proportional to the number of rings in the tube, which is inversely 
proportional to the pitch. In this case, the slope of the friction factor 
versus pitch should be -1. That this appears to be the case is also 
shown in Figure 2. 
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As the spacing becomes very large compared to the reattachment length, 
two other effects will begin to control the pressure drop. The first of 
these is that the wall friction becomes a significant portion of the 
pressure drop. The other is that rotation of the core flow becomes 
significant. The friction factor may, therefore, be expected to asymptotically 
approach the value for a smooth tube at large pitch. It should be observed, 
however, that the wall boundary layer downstream of an annular ring 
is an inlet boundary layer with friction factor higher than that for a 
fully established flow. The "entrance" effects would not be expected 
to be negligible for 25 or more diameters downstream. 


2) Fully Established Flow Through Tubes Containing Helical Wires 
With p/D Greater than about 0.8 

In a long tube containing only a helical wire insert, some rotation 
of the core flow may be expected to occur, since a component of the 
drag force exerted upon the fluid by the wire is in the direction to 
cause rotation. A simple model of this flow has been developed. The 
model studied consists of a fluid flowing in a tube containing a helical 
wire insert. The helical wire insert causes a torque on the fluid which 
rotates it in the direction of the coil, and skin friction opposes the 
rotation. Making the idealization that the fluid bulk rotates as a solid 
body, but slips with respect to the wire insert, an equation may be 
written for the balance of forces in the direction of rotation normal 
to the tube axis . The wire geometry and force system directions are 
shown in Figures 3a and 3b. 


The force exerted on the fluid by the wire, measured in the direction 
normal to the wire, is 


F 

n 


C_ d SL —Pv ‘ 
D w w 2 n 


(3) 
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NOTE: THE ANGLE BETWEEN 
WIRE AND TUBE AXIS HAS 
BEEN INDICATED IN EACH 
DIAGRAM BY a 


Figure 3a. One pitch length of helical wire, unrolled. 



Figure 3b. 


Force system exerted .on fluid by wire. 



Figure 3c. Velocity components. 


where d , wire diameter 
w 


i , wire length 
w 

C D , drag coefficient 
w 

p , fluid density 

v , velocity component normal to wire. 


i , the actual wire length, is given by 


i - i J i ♦ ft® ) 2 

w a V \p / 


where £ , axial tube length 
a ° 


p, wire helix pitch 
D, tube inside diameter 


The component of F in the direction normal to the tube axis is 

n 


F 



The frictional force in the direction normal to the tube axis is 




it B i 

a 


( 6 ) 


where £, friction factor (0. 02 used in the calculations) 

v^ , velocity component normal to tube axis, at inside tube wall, 
neglecting boundary layer effects 

Where the flow is fully established, so that v t is not changing, the following 
force balance must hold: 


15 



( 7 ) 


F ^ = F 
ft wt 


therefore 


C d v 2 = f 7T D v 2 
D w w n 4 t 


( 8 ) 


or by rearrangement. 


n / J 7 r D 


, 1/2 


4 C_ d 
, D w, 


(9) 


From Figure 3c, it may be determined that 

7 r D 


v = 
n 


P a 


1 + 


(x)' 


(10) 


This relationship may be used t o eliminate v from Equation ( 9 ) to yield 

n 


7T D 
~ 


i + 


^,2 1 ) 1/2 

-) J S * 1 


( 11 ) 


where is mean axial velocity 


k = 


7T f D 


4 C, d 
D w 
m 


'D , a modified wire frontal drag coefficient , is defined by 
m 


C D " X * 41 \ D 
m 


w 


2/7 


w 


( 12 ) 


This modified drag coefficient includes the decrease in velocity of the flow 

impinging on the wire due to the tube boundary layer. The form is based on 

a l/7 power velocity profile. has been taken equal to unity. 

w 
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If the flow did not slip with respect to the wire, the velocity ratio 
would be 



no slip 

The ratio of rotational velocity with slip to that without slip is there- 


slip 


no slip < k 


1 + 


f7T D 


2] )l/2 


(14) 


+ 1 


A graph of this relationship as a function of coil pitch and coil wire diameter 
ratio is shown in Figure 4 . 

The axial pressure drop consists of two terms. One of these is the 
skin friction increase due to the angular motion. Neglecting both the change 
in equivalent diameter caused by the presence of the wire and change in fric- 
tion factor due to higher helical flow Reynolds number, the friction term is 


A?, 


„ a 1 2 
£ ~ ~ p v 
D 2 a 



(15) 


where £,, helical flow path length along surface 
h 

v , surface helical velocity with slip (neglecting boundary layer effects) 



a 1 2 

The quantity £ — — p is of course the pressure drop for an axial 

flow through a smooth tube. 
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1 1 1 

1 1 1 1 1 1 1 1 
j ^- = 0.143 

1 


^^= 0.0634 





The second term in the axial pressure drop is due to the axial compo- 
nent of drag on the helical wire. Normalized on the same basis as the friction 


term, this component is 


drag 


(t) 


A? 


smooth 


L[i ♦ (^f 

i i/2 i l 

1 l + 1 

< L V p / 

J ) 


(16) 


From geometry, 


i v 

_h = _h 

i v 

a a 




1/2 


(17) 


The pressure drop in a tube containing a long pitch helical wire may be written 


^ P h w ^ 

[ 

l 1 + 

\ 

m 2 i 

3/2 

(t) 

3 

A) 

smooth j 

L ( 1 + (-)T + 1 

i 7 

Jl + (— ) 2 ) 

- 1/2 



L \ ' p ' )\ 1 

1 ( 

\ v p / / 



A graph of this function, using values of k appropriate to helical wires for 
which friction factor data are available is shown in Figure 5. 

It may be seen that the analytical curves have the same shape as the 
experimental data, and that the proper relationship between wire diameter 
ratio and friction pressure drop is predicted by the analysis. For the larger 
diameter wires the agreement between analysis and experiment is very good. 
The pressure drop for small diameter wires is overestimated by the analysis. 

As can be observed in Figure 4, the slip becomes quite large for small 
diameter wires and short pitches , and it might be expected that departure from 
solid body rotation, disturbance to the velocity profile, (which probably bears 
little similarity to a l/7 power velocity distribution at short pitches) and in- 
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Figure 5 . Comparison of helical wire 
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creased amount of the surface area involved in the recirculating stalled areas 
downstream of the wires, all contribute to cause departure of the actual flow 
from that modeled. One would expect also that flows through thick wire coils 
which allow less slip at a given pitch than thin wires would more closely approx- 
imate the model flow, and this does appear to be the case. 

The analysis should be adequate for making pressure drop calculations 

for wires of d /D > 0. 06 at p/D > .8. Additional discussion of the compari- 
w 

son of Sams’ data and the analysis may be found in the next section. 


3) Transition From A Wire-Wrapped Plug to a Helical Wire Insert 

Consider the single phase fluid which flows through a tube with an abrupt 

change in insert geometry. The insert geometry changes from a helical vane 

to a helical wire coil at the station S. =0. The pitches of the two inserts need 

a 

not be the same. The asymptotic case to which this transition flow tends far 
downstream is given by the preceeding section of this report. The initial 
rotation is known since the fluid is constrained to follow the helical vane, 
and resembles solid body rotation. Thus, the transition flow may be modeled 
as a flow with solid body rotation which decays from one known value to another 
known value. The damping forces causing the rotation to decay may be 
extracted from the asymptotic solution represented by the helical coil analysis 
of the preceeding section. 

Consider the rotational accelerations on a mass of fluid rotating as a 
solid body. The equation of motion is 


M = I 


dw 

dt 


(19 ) 
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where M, applied torque 

I, rotational moment of inertia 
w, angular velocity 
t, time 

Since the axial velocity is a constant through the transition process, 

axial distance and time are linearly related, so that 

dJL 

dt = — (20) 

a 

where £ , axial length 

ft 

v , axial mean velocity 

ft 

The angular velocity may be related to the velocity component paral- 
lel to the surface and normal to the tube axis so that 

2v e 

w = — (21) 

where v , velocity component normal to the tube axis 

o 

D, tube diameter 

The moment of inertia of a fluid mass one pitch length long is 

I = p p - D - (22) 

32 

where p, fluid density 

p, helical wire pitch length 

The torque applied to the fluid mass is equal to the difference between the wire 
drag force component exerted in the direction of rotation, and the wall shear 
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force component exerted opposite to the direction of rotation, both multiplied 
by the inside radius of the tube. From Equations (3), (4), (5), and (6) (proceeding 
Section) this net torque in the direction normal to the tube axis is equal to 


M = ( C d p ~ p v 2 
V D w * 2 H n 


£ 2 
g P * Dp v t 


From Equation (10) 

(24) 

Equations (19), (20), (21), (22), (23), and (24) may be combined to yield a 
differential equation relating v and £ . 

t cL 

This equation is 

2 Hy 

P P l P v a _L (25) 
8 d£ 

a 

Define the following normalized variables: 

£ 

I =-?- 

P 




V t 



f 

where v t , asymptotic value, given by Equation (11). 
f 

By use of Equation (11) , 



where k , defined as before a 


7T £ D 



(26) 


23 



Introduction of the dimensionless variables and rearrangement yields 



8C d 

D w 


1 + 


(?)] 



( 27 ) 


This differential equation is of the form 


dy 


a V. 


+ bT f + 


c = - 


t 

d£ 


( 28 ) 


The equation is separable and easily integrated. If l is measured so that the 
insert type changes at F = 0, then the initial condition is 

F = 0, V = V t ( 29 ) 
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After integration and considerable algebra, the result may be written in 
the form 




For the special cases in which the helical vane and the helical wire inserts 
have the same pitch. 




1 


(31) 


= M + 1 (32) 

This result is taken directly from Equation 14. The solution 
then reduces to 

-k 

(1 + A*) - (1 ~ /£) e 

-k 

(1 + H) + (1 - H) e 

A graph of this function is shown in Figure 6. 

The characteristic length for this transition is 

£ = Y* (34) 

Two graphs of this function are given here. In Figure 7, j^/p is shown as a 
function of wire diameter and helix pitch. In Figure 8, ^zf/D is shown. For 
the longer pitch helices, fewer turns are necessary for the flow to accommo- 
date to the new conditions, but a greater tube length is required because of 

the lengthened pitch. 

It is interesting that many diameters of tube are necessary for flow 
accommodation. It appears that for long pitches and small wire diameters, 
fully established flow would be rather difficult to attain. 

The tube which Sams (Reference 7) used to obtain his data for helical 
wire systems was 0.41 inch ID and 24 inches long for a t/D ratio of 59. This 
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Figure 7. Characteristic length as a function of wire diameter and 
pitch, normalized using pitch length. 









length is barely adequate to obtain fully established turbulent flow in a smooth 
tube (Reference 11) and for the longer pitches he tested was probably inade- 
quate in length for fully established rotational flow. Thus, one might expect 
his data to correspond to cases in which the flow is changing along the length. 
Unfortunately, he does not publish axial temperature profiles in which a vari- 
ation, if it occurred, might be seen. It is certainly true that checks he made 
with no insert agree with published data for fully established flows. The in- 
let to the test section in his apparatus appears to have been a sharp-edged 
reduction, which because of the very high turbulence level generated, would 
tend to reduce the entrance length effects and make the smooth tube results 
agree more closely with published data, but this high inlet turbulence level 
would probably not enhance the rotational component accommodation signifi- 
cantly. 

In the discussion of Section A2, the discrepancy found between pre- 
dicted helical wire friction factors and Sams* experimental data for small 
wires was attributed to departure of the actual flow from that of the model. 

It is likely, however, that a portion of this discrepancy is actually due to the 
fact that flow was not fully established in the experimental apparatus. 


4) Pressure Drop in the Transition to a Helical Wire Insert 

The helical vane to helical wire transition results obtained in the 
preceeding section may be used to compute the integral average friction factor 
for a single phase flow through the transition. This is accomplished as 
follows . 

First the equations for helical wire friction pressure drop are re- 
written in terms of the rotational velocity component v^. After transformation, 
these equations are as follows. Increase in friction due to longer helical path: 



ia 1 2 

D 2 P 



(35) 
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Pressure drop due to drag on helical wire 


4?, 


»dA 


1 2 
f "^ V a 


(36) 


The sum of these two terms, which is the normalized axial pressure drop, 
may be used to define a friction factor for helical wire flow, so that 




nw 




axial 


’hw 
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(37) 


= K 

In the transition flow v^. is a function of axial length from the transi- 
tion. If it is assumed that the local friction in a changing flow is no differ- 
ent from that at the same local conditions in a fully established flow (a com- 
mon assumption for varying flows) then a mean value for the helical wire 
friction factor may be defined. 
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^av ~ 1 f l 
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The velocity component v^ in the transition from a helical vane to a helical wire 
of the same pitch is given by Equation (33) . By definition 
V* 


'f 


(40) 
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v v. 
a 


(41) 


'i 

The second ratio is a constant for a given wire system and is given by 
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( 42 ) 
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With the relationships computed above, the average friction factor, which is 
a function only of the geometry and tube length, may be computed for a given 
case. Evaluation of the integrals must be done numerically, but may be ac- 
complished with a computer. Once the average friction factor for a given 
geometry is computed, the helical wire transition section may be treated as 
a section of smooth tube. 

5) Heat Transfer in the Transition to a Helical Wire Insert 

The single phase heat transfer conductance for helical flow is given 
by Equations (216), (217), and (218) in Section E3. These relationships are based on 
the idealization that the forced convection effects and the free convection due 
to radial acceleration may be added linearly. For a helical wire system all 
of the equivalent diameters may be taken as equal, as a reasonable approxi- 
mation. /?, the volumetric thermal expansion coefficient is equal to l/T for 
a perfect gas, where T is the absolute temperature. Re^ is the linear flow 
Reynolds modulus, given by 


Re * 


NUj. 


inear turbulent vapor 


v D p 
a e v 



(43) 


= 0.023 (Re £ ) 0,8 (Pr) 0,4 (44) 

From these relationships, and the concepts employed in the computation of 
the average friction factor in a varying flow, it may be shown that the average 
Nusselt modulus is given approximately by 

^helical wire 

Nu 

, linear turbulent flow 



( 45 ) 
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where B = 


7.1 (/? AT) 1/3 
0. 13 0.07 

R6 i Pr 


Again, the integrations must be performed numerically, but they depend only 
on the geometry of the insert and its length, and are independent of the flow 
properties as long as the flow is single phase and turbulent. 
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B. The Motion of Droplets Entrained in a Vapor Flow 


In the general case of forced convection boiling with entrained liquid, 
droplets of time- and space-varying size accelerate in a vapor stream which 
is non-uniform in three directions. A gravitational field, rotation about the 
tube axis, possible collisions between droplets or between droplet and wall, 
and droplet number density change by ejection from a liquid layer, may all 
add to the complexity of the motion in a general analysis. 

It is clear that a single comprehensive treatment of droplet motion, even if 
suitable functional relationships were known, would be very difficult to interpret, 
and would be far too inclusive for many practical problems. In order to maintain 
the maximum utility in these analyses of this section, many restricted models of 
various aspects of entrained droplet motion have been analysed. Several models 
are also presented which are developed by synthesis of the simple models to 
study a complex motion. 

1) Droplet Number Densities in a Boiler Tube Containing a Fog Flow 

In a boiler tube with a dry wall at qualities less than unity, the number 
density of the droplets entrained in the vapor depends upon vapor quality X, 
vapor density p liquid density p vapor velocity v^, droplet velocity v , 

droplet diameter D , and tube diameter D. If the droplets are idealized as 
P 

spherical, the relationship between the variables may be shown to be 

N = 1.5 _^v 1 (46) 

\ p i (v /v ) (D /D)3 

p v p 

where N is the number of droplets in an axial length equal to the tube diameter. 
This expression has been evaluated and is presented in Figure 9. 
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Figure 9. Number of droplets in an axial length equal to one 
diameter. 
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2) Axial Velocity Component of Constant Mass Droplets in a Uniform 


Velocity Stream 

This analysis pertains to the motion of entrained droplets in an upward 
flowing vapor stream. In order to eliminate heat flux as a parameter in the 
analysis, the idealization is made that the change in vapor mass flow rate is 
small over the boiler length of interest. 

Other idealizations are: 

1) droplets are of constant mass 

2) droplets are spherical 

3) the flow is vertically upward so the gravitational field opposes 
the motion 

4) the droplet drag coefficient is constant* 

5) the liquid phase occupies negligible volume in the boiler 

6) the pressure gradient is negligible 


*For commentary on this simplification, refer to pages 52 and 55. 
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7) the axial component of vapor velocity is uniform across the boiler. 


The normalized equation of motion is 


dV P 

dt 

a 


d - V 2 - % 


where 


v 

P 


v 

J? 

V 

a 


t 

a 


t 


T 

a 


or r 

- e a 
~~ ”v 

a 

v = droplet velocity 
v^ = vapor axial velocity 

t = time, measured from initiation of droplet motion 
g^ = earth 1 s gravitational field 


r v f. atlmescale 

a \4 Dp „ t a / 

= droplet drag coefficient, numerically equal to 0.4 

D = droplet diameter 
P 

= vapor density 
Pp = liquid density 

This equation may be integrated readily and rearranged to yield the relation- 
ship between velocity and time 


(47) 
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V p = (l-g e >- 


1 -e 


(1 - NTgJ - (1 


2 \[g~ t 
&e a 

g)e 

The relationship between axial location and time is obtained from a second 
integration. 

_ dz 
V P = dt 

a 

where 


(48) 


z = 


V T 

a a 

z is the axial length from initiation of droplet motion. 

By integration of Equation (48) and rearrangement one obtains 

2 \[g t 

(w iA 


y e a 


t 1 + N^e)' 


t - lo<; 
a 




2n/i: 


(49) 


Under many circumstances, where the axial drag force is large and small axial 
distances are of interest (so that droplet slip velocities are much larger than 
the asymptotic slip due to gravity) the gravitational force term is negligible. 
The resulting simpler form of Equation (47) , namely 


d v 

P _ 2 

-Tf- - (1 ‘ V ) 
ci t» a 


(50) 


may be integrated twice to yield 


and 


v = 1 - = 

P 1 + t 


z = t a - log (1 + t a ) 


(51) 


(52) 


The results of Equations (48), (49), (51), and (52) are shown in Figures 10 and 11. 
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Figure 10. Droplet accelerating in a uniform stream: Axial velocity 
versus time- 
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3) Axial Velocity Component of Constant Mass Droplets in a Boiler 
Tube Subjected to a Uniform Heat Flux 

This analysis is of a model of droplet motion in a vertical upward 
flow in which the vapor flow increase due to vaporization is significant within 
the boiler length of interest. 

The idealizations on which the model is based are: 

1) droplets are of constant mass 

2) droplets are spherical 

3) the boiler is subjected to a uniform heat flux 

4) the flow is vertically upward, so that the gravitational field op- 
poses the motion 

5) the droplet drag coefficient is constant 

6) the liquid phase occupies negligible volume in the boiler 

7) the droplets are injected with zero axial velocity at an arbitrary 
location in the boiler 

8) the vapor velocity is uniform across the tube diameter. 

These idealizations lead to the following normalized differential equation 


dv 




r - .2 _ p 

0 - v \ - y = v — - 
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p d X 


where 
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Se 
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A X 
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X 
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a characteristic length 


v 


r 


A X^, a characteristic velocity 


(53) 
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where 


X = * reciprocal of a characteristic time 


= drag coefficient of droplet, taken equal to 0.4 
D = tube diameter 


D = droplet diameter 
P 

g e = acceleration of gravity of earth 
£ = axial length 

£ q = point of injection of droplet 
L = latent heat of vaporization 


— = heat flux 
A 


v v = vapor velocity 
^ = droplet velocity 

y = normalized gravitational force 
In this coordinate system the vapor velocity at any point is given by 


v 

v 



= £ 


r r 

Equation (1) has been numerically integrated for several droplet injection points 
£ , and for several values of the gravitational parameter y . The results are 
shown in Figures 12 and 13 . 

The normalized velocity difference (between the vapor velocity and the 

particle) for a particle with suitably long path length in the boiler, defined by 

v - v 


V 

s 



(54) 
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£ 


- v 


P 


£ 


42 







continually decrease^ The variation of this quantity is shown in Figures 14 and 15. 

The variation in spatial location and particle velocity with time must 
be extracted from the above results numerically, by the relation 



(55) 


where 



t = time from instant of droplet injection 
and A = denotes a small increment. 

The initial variation of £ with time is shown in Figure 16. 

Because of the gravitational field, there is a minimum vapor velocity 
below which upward motion of droplets will not occur. This minimum is given 
by the relationship 

J 2 = y (56) 

o 

By use of the definitions of the normalized variables, Equation (56) 
may be reduced to 


a 

o 


1 L 



/Vv gD p 

V C D 


(57) 


Droplet Velocities in a Specific Boiler 

The boiler chosen for this example is described by the following 
parameters . 

Fluid: Potassium at 1800 °F 
Tube Diameter: 0.42 inch 

2 

Heat Flux (Uniform): 28 Btu/ft sec 

2 

Mass Velocity G: 30 lb/ ft sec 

These conditions correspond to some experimental conditions reported in Reference 2. 
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Figure 14. 
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Constant mass droplets accelerating in a uniform heat flux boiler: Minimum (asymptotic) 

slip versus location in boiler. 






The particular point in this boiler to be considered is 
Quality X: 0.8 

■Length from Onset of Boiling: 5.73 ft. 
Computed parameters of interest are 

Jf) 

A p D L 
v 

27.7 

sec 
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Because actually observed droplets may not have been in the stream long enough 
to attain minimum slip, one would expect droplets of these. sizes to have greater 
slip than that tabulated, but not less. 

It may be seen from Figure 15 of this report that the effect of the gra- 
vity term is small for all of these cases. For this specific boiler, entrainment 
of even the largest droplets considered would be expected for locations in the 
boiler above 


JL 

o 


D 


6.9 


This is approximately three inches after the onset of boiling. 


4) Axial Velocity Component of Fragmenting Droplets in an 
Accelerating Stream 

It has been shown in a number of investigations (see for 
example, Reference 12) that a liquid droplet is unstable when its Weber number 
becomes greater than a fixed value, usually taken to be about 16*. The un- 
stable droplet fragments very rapidly, so that the maximum droplet size under 
a given set of conditions may be near that given by the critical Weber number. 
Consider a droplet of varying size accelerating in the vapor stream of a uniform 


*It is recognized that there are no doubt other mechanisms that can play a role 
in controlling droplet size, such as spallation or fission when droplets contact the 
boiler tube surface. In the present analysis an attempt is being made to investi- 
gate drop slip velocities in the bulk region on the premise that the Weber criterion 
is in operation. 
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heat flux boiler. At the instant before fragmentation, the droplet is acceler- 
ating under the influence of the drag force exerted by the accelerating vapor 
stream. After the (instantaneous) fragmentation process, each of the frag- 
ments continues acceleration under influence of the vapor drag- force, but now 
there is a new frontal area and mass for each fragment. On the average the 
fragmentation process does not contribute to the axial momentum of any of 
the fragments. This process may be idealized as continuous. The equation 
of motion is 

*=* <68) 


when 

F is the force applied to the system 

m is the droplet mass 

v is absolute axial velocity of the droplet 
P 

t is time, 

— is the mass change due to fragmentation 
dt 


dxxx 

The term(v^+ Av^-^~~ is the momentum carried away by the fragments. Since 

there is no net axial velocity change due to the fragmentation process, Av^ = 

When this equation is expanded, the two terms containing -fo - ft- , which have 

dt 

opposite signs, cancel. The equation is therefore 
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(59) 


F = m 


dt 


For the present purpose, distance is a more suitable coordinate than is time. 

dv dv , . 

0 d cd 

dt ~ dl &t 


dv 


= v 


p dt 

The equation of motion may be written 

dv 

_ P 

F - “vzr 


(60) 


where Si is the axial distance measured from onset of boiling. The drag force is 
given by 

Ok -i r\ 

(61) 


„ „ 7r 2 1 2 

F_ as C_ D — -p V 

D D 4 p 2 v s 


where 


C_ is the drag coefficient 
D 

D is the particle diameter (the particles are idealized 
P 

as spherical) 

p is the vapor density 
v 

v g is the slip velocity. 


v is defined by 
s 


v = v - v 
s v p 


(62) 


where is the absolute particle velocity, v^ is the vapor velocity. The force on 

the droplet due to gravity will be neglected in this analysis. 

In general, the drag coefficient is a function of the droplet Reynolds 

number, but as will be shown later, the Reynolds numbers for droplets of inter- 

3 5 

est in this investigation are in the range of . 5 x 10 to 2 x 10 , where the 
drag coefficient is almost constant at = 0.4. For the present purposes it 
will be assumed that the drag coefficient is invariant. Neglecting the liquid 
volume, and considering the case in which the boiler is subjected to a uniform 
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heat flux, the mean vapor velocity at a cross section is linearly dependent 


upon length f . 


where 


v» 

v 




p D L 
v 


-r- is the heat flux 
p^ is vapor density 
D is tube inside diameter 


( 63 ) 


and 

L is latent heat of vaporization. 

On an average basis one may calculate the droplet diameter from the 
Weber number, so that 

D 


16<7g 


p v 
v s 


where 


cr is surface tension of the liquid 

D is droplet diameter 
P 


The droplet mass may be calculated from 

7T 3 

m =-~p„ D- 
6 JL p 


where 


(64) 


(65) 


PjL is the liquid density. 

By combining Equations 60 through 65 and rearranging, one obtains 


3 C D p v 2 = V p dV p 

64 crg c PjL V e - 


df 


( 66 ) 
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The slip velocity may be replaced by 

v = v - v 
s v p 

- 4f) 

p D L 
v 

Define the reference length and velocity 


Z-v 


4 tr Vi D 




V =' _ _ - 

r p D .L r 
v 

Also define the dimensionless variables: 
partiole velooity 


V = — 1 
P ▼ 


axial length 


s =• 



_ 1 _ 
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(67) 

( 68 ) 
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In this coordinate system, the normalized vapor velooity, defined by 


V = — 
v v 


has a slope of one. Equation (66) becomes 


dV 
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<b > V p ) 


4 ds 


(72) 


or 


dV 
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This equation is nonlinear and apparently cannot be solved in closed form. 

An approximate solution has been obtained by the method of isoclines 
for several values of starting position for the droplet. These results are shown 
in Figures 17 and 18 . It may be seen that the droplets accelerate very rapidly 
and approach a limiting slip. Define a normalized local slip velocity as before (Eq. 54) 

v v - v 

g. = V P (74) 

V V - 

V V 

= V 
s 

The local slip velocity for a particle which has been in the stream long enough 
to^come to equilibrium acceleration varies as shown in Figure 19. According 
to this analysis, the droplets fragment and accelerate so rapidly that one would 
expect most droplets entrained in a vapor flow at a particular axial station to 
be traveling at nearly the same velocity regardless of whether they had been 
entrained early in the boiler or only a short distance upstream from the point 
of observation. 


5) Size of Fragmenting Droplets in an Accelerating Stream 


In the calculation of slip velocity with acceleration, it has been as- 
sumed that the droplet size is defined by the maximum Weber number for 
which liquid droplets are stable, namely Equation (64). Equation (64) may 
be transformed by use of Equations (63) and (74) to yield 



(75) 


Introducing the characteristic length defined by Equation (69) one obtains 
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This relationship is shown graphed in Figure 20 . It may be seen that the de- 
crease in droplet size is very rapid during the initial portion of the boiler. Be- 
cause of the nature of this analysis, the predicted droplet size is near the max- 
imum one would expect to observe. An unstable droplet probably fragments in- 
to several particles all of which would be below the critical Weber number. 
However, if a droplet fragments into for example two particles of equal mass, 
the diameters of these two particles are 80% of the diameter of the original 
drop, or almost the same size. Therefore, unless fragmentation is into a 
very large number of droplets , the observed diameter range would not be ex- 
pected to be very great due to this effect. 

The velocity of the vapor is not uniform across the tube, and not all 
droplets are subjected to the same accelerations. Although accelerating flows 
have more uniform velocity profiles than fully established flows, the liquid 
vaporization on the wall injects vapor into the core flow with no axial velocity 
component, and this would be expected to cause a velocity profile distortion 
near the wall. Therefore a scattering of droplet size due to the nonuniform 
velocity profile would be expected. 

The idealization was made that the sphere drag coefficient is constant 
for the range of droplet Reynolds numbers of interest here. The Reynolds num- 
ber, based on the slip velocity, is 


Re 



(77) 


This relationship may be transformed to 



For Freon 12, typical conditions are (see Reference 12) 
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Figure 18. Normalized velocity of accelerating fragmenting droplets, s greater than unity. 
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The characteristic length is l = 0. 55 ft and the parameter in the Reynolds 
number has the value 6. 9 x 10 . The Reynolds number range is 


s 

1 

10 

30 


Re 

c 

2.4 x 10 c 

r 

. 93 x 10 
. 66 X 10* 3 


These Reynolds numbers are within the range for which is nearly constant. 

The results of this analysis predict very much different behavior for 
a low surface tension fluid such as Freon, and a high surface tension fluid 
such as mercury. Typical data taken from actual experiments are 

M ercury . (I0Q°.Q .Freon. Ulfil EL gfitaasiam I3MIB 

Heat Flux , 2 BtU ' 2 10 28 

ft sec 

Tube Diameter, inches 0.18 0.452 .42 

These values lead to the following computed conditions at a position 8 feet from the 
initiation of boiling: 


Droplet TSize Parameter 

0.096 

.208 

.16 

(in Equation (76)) 




Reference Length ft 

12.6 

.548 

4.6 

V at 8 ft 
s 

0.77 

.13 

.6 

Vapor Velocity , ft/ sec 

114 

41.9 

216 

Droplet Velocity , ft/ sec 

26 

36.5 

108 

Droplet Diameter, inches 

0.057 

.027 

.08 
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6) Deposition of Entrained Droplets in Linear Flow 

The vapor flow in a linear flow boiler tube almost always has 
some liquid droplets entrained in it when the quality is less than unity. If the 
wall temperature is not so high that film boiling occurs, and if the liquid wets 
the boiler tube wall, then presumably liquid droplets which strike the wall spread 
out and vaporize. If there is a continuous liquid film on the wall, the impinging 
droplet would be expected to become a part of the film and not to be returned to 
the vapor flow. It is possible of course, that entrainment of liquid from the film 
may also be occurring, but the idealization will be made that these processes 
are independent and may be treated separately. 

The rate of deposition may be expected to be related to the 
concentration of droplets entrained in the vapor, and if the droplets are suitably 
small it may be valid to treat the deposition process as a process of diffusion 
of very small particles from the bulk to the wall. The controlling equation is 
then of the form 


w 

where d 

A 



c 

s 



h (c - c ) 
d s w 


is the mass flux impinging upon the tube wall 
is a diffusion conductance 


is the vapor stream droplet concentration 


(79) 


and c is the concentration of droplets at the wall, 

w 

According to the present flow model, the wall acts as a "sink" for the droplets, 

so that c =0. According to Reynolds' analogy between the turbulent transfer 
w 

of mass and heat, (see for example Reference 13) 


h 


c 


P 
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where h is the heat transfer conductance and c is the vapor specific heat. Con- 

P 

tinuing the analogy, the Dittus - Boelter equation for the Nusselt modulus des- 
cribing turbulent heat transfer in tubes is transformed to 


where 


D 

2 > 

Re 


h d° 

Z) 


= 0.023 Re 0,8 ( Sc j** 4 


is tube inside diameter 

is the molecular diffusion coefficient 
is the Reynolds modulus, defined by 


Re = 


V V D Py 


(SO) 


V 

V 


V 

Sc 


is mean vapor velocity 
is vapor density 

is vapor viscosity 

is the Schmidt modulus , defined by 

M.. 


Sc = 




In the present analysis the Schmidt modulus is taken equal to unity. The mean 
liquid concentration in the vapor core is 


where 


w_ 


w e 

c = 

s' v A„ 

P T 

is the entrained liquid mass flow rate 
is the droplet mean velocity 
is the tube cross section area. 


( 81 ) 
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The decrease in effective diameter of the tube due to a possible liquid 

film on the wall has been neglected. The mass flux to the wall, w , may be re- 

d 

lated by continuity to the loss in entrained flow. 


dw 


w 


IT D 


(82) 


di A 

where £ is the axial length, measured from a datum location as yet unspecified. 
Combining Equations (79), (81), and (82), one obtains 


d w 4h , w 

e d e 


df 


v D 
P 


(83) 


In the analysis of Section B.4, a relationship between the droplet velocity v^ and the 
axial coordinate £ was derived for fragmenting droplets in a boiler with uniform heat flux. 
This relationship is of the form 


v = v d> (a) 
p r 


(84) 


where 


£ 



p D L 

v 


c t is liquid surface tension 

C D is droplet drag coefficient, equal to 0.4 for the range of drop- 
let slip Reynolds moduli of interest 
p is liquid density 

L is latent heat of vaporization 

— is the uniform heat flux 

A 

v is the vapor velocity, a function of length, 

v 
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For the uniform heat flux case, and neglecting liquid volume in order to be con- 
sistent with preceding approximations, 


_q 

4 A 

T v " Lp D 
v 


= XI 

then 

dw 4 w i / h , (s)\ 

_ e r / d - ' \ 

ds “ D \v r <p (s ) ) 

the diffusion conductance is given by (80), (85), and (86) and is 


( 86 ) 
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where 
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Combining Equations (83) through (88) , one obtains an equation for the entrained 
flow 


where 


dw 

e 

ds 


k l (f>( a) 

4g 
D X 


(89) 


For cases in which entrainment of liquid into the vapor stream from the wall 
film is negligible, then no liquid is added to the vapor flow along the boiler. For 
these cases the entrained flow has some point at which the liquid droplet flow is a 
maximum. This will be taken as the initial condition for the entrained flow. The initial 


condition is therefore 


(89a) 



Integration of Equation (89) subject to (89 a ), yields 


w 


w 


= e 


K ---1 ,, 

A = e 1 


( 90 ) 


The dimensionless constant in the exponent is 


LP 


k = .0758 — 

1 q D 


j L / Jr 2 \ 3 

L C D / 


0.8 


(91) 


It remains to evaluate the integral I. 


The functional relationship between the particle velocity v and s, denoted sym- 
bolically by 4 >, cannot be obtained in explicit form. The slip velocity, defined 
by v g = v y - v , is fit approximately, for s > 1, by the relationship 


v 

s 



(92) 


In Figure 21, Equation (92) is compared to the particle velocity given in Section B.4. 
It may be seen that for present purposes the fit is acceptable. 


v v v 

p _ _v _s 

v ~ v v 

r r r 



■ = *<s) 

The integral exponent in Equation (90) becomes 


I 



1.5 

s 


ds 


s 1 ’-0.87 


(93) 
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in Section B.4. 



For values of s > 1, the integrand is fit acceptably by the expression 


1 1 5.7 

0.2 + 2 + 7 

s s s 


(94) 


A comparison of the integrand in Equation (93) and the approximate relation of 
Equation ( 94 ) is shown in Figure 22. From the integral with the approximate 
integrand of Equation (94), the result is 

t „ , 0-8 0.8 1 1 „ „„ 1 1 

I = (1.25) (s -* 1 )-(--—)- .95 (— - — ) (95) 

1 s s. 


Examples of the predicted droplet deposition rates based on this analy- 
sis have been prepared for two representative cases, one Freon 12, and the 
other potassium. The conditions are 

Freon 12 Potassium 


Specified Conditions 


= 110°F 
= .452" 

Btu 

= 10 — 

*,2 

ft -sec 

Xats. = . 516 
Computed Constants 

S 1 

l 


= 5 
= .548 


T 

D 

q 


T = 1800°F 
s 

D = . 42 " 


_q 

A 


28 


Btu 

2 

ft -sec 


X at = . 488 



1 
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The resulting deposition rate graphs are shown in Figures 23 and 24 . Although 
these curves are not exponentials in the coordinate s, an approximate effective 
characteristic length can be obtained from the variation of 


w 


w 


as a function of s. For potassium, at 

w 


— - “ . (■“■j) = .8 


this corresponds to an axial tube length of 3.68 feet. For Freon 12, at 

W e 1 

W -T’ <S ‘ S 1> - - 8 

e o 


This corresponds tQ an axial tube length of .438 feet. The rate of decrease of 
entrained flow is almost an order of magnitude higher for Freon 12 than for 
potassium at the specified conditions. 

7) Axial Path Lengths for Droplets with Transverse Velocity Components 
Several of the preceeding analyses have described the axial motion of 
droplets in a boiler tube. There are a number of mechanisms by which droplets 
might obtain significant transverse velocities as well. One such possible 
mechanism, which is described and analyzed in Section C.5 of this report, is 
due to bursting of vapor bubbles in a nucleate boiling liquid film. 
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The axial path length of droplets with transverse velocities is studied in 
the following paragraphs, with calculations based on velocities obtained from the 
ejection mechanism. Although the calculations of the mechanism are based on 
water observations, there is a general similarity between the surface tension 
and density of water and potassium, so one might expect ejection velocities of 
the two fluids to be similar. 

Droplets ejected from an annular liquid film in forced convection nucleate 
boiling have a significant velocity component normal to the film surface. If 
the flow in the boiler tube is vertically upward, the gravitational field enters 
into the axial motion only, drag on the vapor is the only force acting in the 


transverse direction, and the transverse equation of motion is 

-o. . dV 


tr 


dt 


(96) 


tr 


where v^is the transverse velocity component and the constant is defined on page 

130. The appropriate initial condition is: t = o; v. = v . This equation is easily 

tr tro n 

integrated to yield the following relationship between velocity and time: 


-C 


1 


t = 



(97) 


A second integration yields 


C i 
1 t 


= 1 


+ v C t 
tro 1 


(98) 


where f is the transverse distance measured from the film surface. 

These relationships may be easily expressed in the notation used for the 
analysis of accelerating droplets in a uniform heat flux boiler (Section B.3) 
Equation (97) becomes 


-t 



(99) 
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where 


and 


t = At 


V v 


Equation (98) becomes 


e = 1 + v. t 
tro 


(100) 


the maximum transverse path length is equal to the tube diameter, so that 


D 

jr = — 

tmax i 

r 


With the specification of ^n initial transverse velocity, Equation (100) defines 
the maximum flight time for a droplet, and the distance the droplets are car- 
ried axially may be taken by interpolation from Figure 16. The results of such 
a calculation, for the same conditions as were used in Section B.3 (page 45), 
is given below. In the absence of any other information on the velocity 
of droplets ejected from nucleate boiling films, the same velocity as observed 
for water droplets, 4 ft/sec was used as a basis for the calculations (see page 131). 


Droplet Diameter, 
inches 

t 

it - je 

o 

A i 
D 

-3 

2 x 10 

.276 

10 

40 

5 x 10” 3 

.255 

3 

30 

-2 

1 x 10 

.249 

1.13 

22.6 

2 x 10 2 

.246 

.35 

14.0 


As one would expect, the smallest droplets are swept the farthest downstream 
before impacting with the opposite wall of the tube. The two mil droplets are 
carried 40 diameters axially. 

The two largest droplets would not have been in the stream long 
enough to have attained the predicted asymptotic slip. In fact, the 20 mil drop- 
let would have only about half the velocity of the vapor at the time it impacted. 
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8 ) The Motion of Large Liquid Droplets in Helical Flow 

It is usually assumed that droplets entrained in a helically flowing 
vapor in a boiler tube concentrate near the wall. However, droplets swept off 
a helical divider may take a significant amount of time to migrate to the wall, 
and since they travel in a spiral path to the wall, slipping with respect to the 
vapor and the helix, they could actually be redeposited on a helical vane rather 
than on the boiler tube wall. During the time they take to migrate to the wall 
the droplets are also swept axially downstream, and the axial flight time is also 
of significance. The study reported here is intended to provide some informa- 
tion about droplet paths in helical flow. 

Idealizations which describe the model upon which the analysis is 
based are listed below: 

1) The droplets are rigid and spherical 

2) The axial vapor velocity is uniform 

3) The rotational velocity field is similar to solid body rotation, 

so that the rotational vapor velocity is a linear function of radius 

4) The droplet drag coefficient is constant 

5) The liquid volume is negligible. 

6) The droplet enters the flow with negligible velocity at a known 

point x , y , z 
o o o 

7) The droplet velocity in the transverse directions is negligible 
compared to the vapor velocity. 

Since the axial vapor velocity is uniform, the motion in the z direction (lying along 
the tube axis) is independent of the x and y (transverse) motions. Therefore, only 
the motion in the x, y plane is considered here. The coordinate system is shown in 
Figure 25. The equations of motion are: 

x. direction 


F 

x 
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( 101 ) 


& 


where 


is the vapor drag force component in the x direction 

m is the droplet mass 
t is time 

x is a transverse coordinate measured from the tube axis. 


y direction 


where 


d 2 y 

F = m — i 
y dt 2 


F^ is the vapor drag force component in the y direction 


( 102 ) 


y is a transverse coordinate orthogonal to the tube axis (z) and x. 
The drag force components are taken to be 


1 2 
F = C A — p (v ) 
x D d 2 v v 


F = C A (v ) 2 

y D d 2 v v 

y 

a 7T _ 2 

A = D 
d 4 p 


(103) 

(104) 


and 


where 


7T 3 

m = ~p D 
6 1 p 


is the droplet drag coefficient 

A J is the droplet frontal area 
a 

is the vapor density 

v is the vapor velocity component in the x direction 
x 


v is the vapor velocity component in the y direction 

y 

0^ is the droplet (or particle) diameter 


Pg is the liquid density. 
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x = r cos 0 
y = r sin0 



Figure 25. Coordinate system. 



Figure 26. Vapor velocity components seen from the -Z direction ( right 

handed helix). Helix and vapor rotation are in the +0 direction 



The vapor velocity at any point is v^ = kr, where r is the radius from the 
tube axis, k is a constant, and the velocity is in the direction normal to the 
radius. 

Consider motion of the particle in the first quadrant, facing upstream 
(in the -z direction) so that the vapor rotation is counterclockwise. The drop- 
let is at point p. A diagram of the vapor velocity components is shown in 
Figure 2. The signs of the drag force components are given in the following 
table: 


Quadrant 

F 

X 

F 

y 

i 

- 

+ 

ii 

- 

- 

m 

+ 

- 

IV 

+ 



The vapor velocity components are given by 

v = -krsinfl 
v 

x 

= -ky 

v = krcosfl 
v 

y 

= kx 

The constant k has the units of (time) , so that the characteristic time r may 
be defined as 

1 


A convenient characteristic velocity is therefore 

a 

T 

V = 


The same characteristic length used in Section B.3 appears here. It is 
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The normalized variables are defined as 



r t 

r 


With these definitions, the equations of motion (101) and (102) become 


d 2 x 


-2 

dt 


r 



(105) 


and 


r 

for the first quadrant. 

These two equations are coupled and non-linear and no obvious method 
of obtaining a closed-form solution is available. They have been solved numerically 
by a tedious but straightforward iteration technique, with a single parameter, x^, the starting 
location. Since the vapor velocity field is axially symmetric, no loss of gen- 
erality results from choosing initial locations along the x axis. A typical trajectory and 
useful nomenclature are presented in Figure 27. The characteristics of one specific 
trajectory are given in Figures 28, 29, and 30. Because of the nature of the dynamic 
model, one might expect that the shape of the particle trajectory in a plane normal to 
the tube axis would be independent of the starting location. That this is the case is 
demonstrated in Figure 31. In this figure, all three normalized trajectories in the 
r, 0 plane have been superimposed and it may be seen that aside from some scatter 
(probably due primarily to inaccuracies in the calculation method) the trajectories have 
a similar shape. The particles from different starting radii do not, however, reach 
similar points along the trajectory at the same times. The variation of the angle# 
with time for the three starting locations is shown in Figure 32. 


_2 

= +x 


(106) 
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Figure 27. Schematic diagram of droplet trajectory in the directions transverse to 
the tube axis. Origin of coordinates fixed at the tube centerline. 

(Projection of helical droplet path onto a plane transverse to the tube axis.) 



Figure 28. Motion of a droplet in helical flow, x 0 - 0.1. 



NORMALIZED RADIUS 



NORMALIZED TIME 


Figure 29. Radius versus time. Droplet in helical flow, x Q = 0.1. 
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, ANGLE IN DEGREES 






0 20 40 60 80 100 120 140 160 180 200 

0 , ANGLE FROM STARTING POINT TO PARTICLE LOCATION 

Figure 31. Comparison of computed particle trajectories for 
three starting locations. 
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9) Droplet Trajectories in a Helical Flow Boiler 

The analysis of transverse motion of entrained droplets in a helical 
flow given in Section B.8 may be combined with any of the analyses of axial motion 
to compute three-dimensional trajectories. As an example of the procedure, in the 
following paragraphs the transverse motion has been combined with the axial motion 
analysis of Section B.3 (constant mass droplets, uniform heat flux). In order to 
illustrate the calculation method, a specific example has been carried through in 
detail. The boiler tube and its operating conditions are identical to those used in 
Section B.3 but it is assumed that a single vane helical insert with a centerbody is 
in the tube. A sketch of the three dimensional motion is shown in Figure 33. The 
conditions are: 

Fluid: Potassium at 1800 °F 

Tube Diameter: 0.42 inch 

, 2 

Heat Flux (Uniform): 28 Btu/ft sec 

, 2 

Mass Velocity G: 30 lb/ ft sec 
p 

Helical Insert: g- = 3.1416 

Quality X: 0.8 

The calculation proceeds as follows: 

vapor flow rate = GX. 

= 24 lb/ ft 2 sec 

axial vapor velocity (uniform) 160 ft/sec 

7T D 

vapor velocity normal to tube axis at the tube surface = v . , 

p axial 

= 160 ft/sec 
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From Section B.8, the relationship for the vapor velocity normal to the tube 
axis is 


where 


= kr 


2 2 2 
= (x + y ) 


Therefore 

_ 160 ft/ sec 
0. 21 inch 

= 9150/sec 

Consider 1 mil diameter droplets. 

From the definition of the characteristic length l , 


i 


r 


D 


40 P. 

p SL 


3 D 



where 

D is the particle (droplet) diameter 
P 

D is the tube diameter 

C D is the droplet drag coefficient (taken as 0.4) 
is the liquid density 
p^ is the vapor density. 

For the conditions stated. 


SL 


r 


D 


= 2 


The droplet for this example is injected into the stream at a radius 

, D 

1/4 of the tube radius (possibly by being stripped from a centerbody of diameter — ). 


88 


Thus 


D 


x = r = „ 
008 


and 


r o / £ 


(*)■ 


1 _ 

16 

The droplet collides with the wall at the point 
r = 4 r 

o 

1 

4 


From Figure 31, 

From Figure 32, by graphical interpolation to r 


= = 4 at© = 118°. 

r 

o 


0.063, t = 11.5 
r 


t = 


t T 
r r 


= 1. 25 milliseconds 


From Section B.3, 


t = At 
a 


* In this section the normalized time t is used for the radial motion. The symbol t 

r a 

is used for the axial motion computed in Section B.3. The two quantities are defined 
using different characteristic lengths. These two time scales must therefore be kept 
distinct. 
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where 


and 


X 


.iill 

p D L 
v 



is the uniform heat flux 


L is the latent heat of vaporization 
For the present conditions, 

X = 27.7/sec 


therefore, 

t = . 0346 
a 

The axial normalized starting point in the boiler is 



(M 


The axial and radial motion characteristic lengths are identical. 

The droplet injection point is at quality X = 0.8, which occurs 140 
tube diameters from the onset of boiling with this heat flux. Therefore 


i =70 
o 

By interpolation from Figure 16 of Section B.3, for the conditions £ = 82 

and t = 0.035 
a 

I -I =1.5 
o 

The axial distance the one mil diameter droplet travels before impact 
with the wall is therefore 


£ - j t £ 

-5-^= »'Vf 


= 3 
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The droplet travels three tube diameters downstream, and rotates 
an angle of 118° about the tube axis. Its distance from the tube axis, which 
started at D /8 has grown until it impacts with the tube wall. 

This droplet may or may not impact on the helical vane. If its starting 
point is in roughly the first third of the distance between two successive pitches 
it will probably hit the wall. If the helix has two vanes (or if there is a multi- 
slot insert) the droplet will almost certainly hit one of the vanes, regardless of 
the axial location of its starting point between pitches of the helix. 

A number of other specific cases of droplet motion have been worked 
out and are presented together with the detailed example in Table 1 . Graphs of the 
effect of droplet diameter on axial length to impact and the effect of helix pitch 
length on axial length to impact have been constructed from this table and are 
presented in Figures 34 and 35. 
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TABLE 1 

With the exception of droplet diameter and helix pitch, all conditions are identical to those 
given on page 86 of this report. 



(1) characteristic time for radial motion 

(2) characteristic length for axial and radial motions 

(3) time to impact with wall (actual impact may be on helix) 


















10) Spacing of Sections of an Interrupted Insert 

In order to remove liquid from the centerbody of a boiler helical 
insert, it has been proposed to interrupt the insert for an interval so that the 
liquid flowing along the centerbody would be shed as droplets and returned to 
the heat transfer wall. The analysis presented here is intended to provide in- 
formation about the proper axial lengths for such insert gaps. 

Several phenomena are taking place simultaneously which affect the 
droplet motion. The droplets, initially at rest, slip with respect to the vapor 
in both the axial and rotational motions, and are accelerated in a three dimen- 
sional spiral path due to drag forces. At the same time the vapor rotation, 
no longer maintained by the helix, begins to decay due to boiler wall friction, 
and the vapor axial velocity continues to increase due to vapor addition from 
vaporization. 

In addition, the droplets may be undergoing fragmentation due to 
internal flow instabilities, vapor shear gradients, or other forces. In order 
to analyse this complex two phase flow, three major restrictions have been 
made, which are: 

1) The axial vapor velocity change is negligible within the axial length 
required for the droplets to impact on the wall. 

2) The vapor angular velocity change is negligible within the axial 
length required for the droplets to impact on the wall. 

3) The droplet size is invariant during flight. 

With these restrictions, the motion of the droplets may be synthesized from 
the helical motion computations given in Sections B.2 and B.8. A sketch of the 
geometry of the composite problem is shown in Figure 36. 
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Figure 36. Sketch of droplet motion. 
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The analysis in Section B.8 of the motion of droplets in a helical 
flow yields the normalized droplet radius from the centerline 
as a function of angle (measured about the centerline from the position droplet 
motion is initiated) and normalized time, with initial normalized radius as the 
single parameter. Because of the nature of the analysis, the helical motion 
utilizes the same reference length as that which appears in the analysis of 
Section B.2 but a different reference time. Denote the radial normalized time by 
t^. The relationship between the two time scales is 

t = t T = t T 
r r a a 

or, by substitution of the characteristic times for the two cases 


t 

r 


P 


2V 7T 
a 


t 

a 


„ C_ p 
1 v 
4 D p 

p & 


v 

a 


t 

a 




(107) 


Droplet helical motions were computed for three different initial radii 
and the trajectories were shown to be geometrically similar in Figure 31 of 
Section B.8. The relationships between angle and time and between normalized 
radius and time, not presented previously, but cross-plotted from Figures 31 
and 32 are shown in Figures 37 and 38. In these figures, the normalized time is 
expressed as 

t 

r 

t 

r l80° 
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where t is the normalized time at which the droplet has traveled halfway 

r i80° 

around the tube from its starting point. Figure 39 shows the variation in time 
to travel halfway around the tube as a function of centerbody radius . This re- 
lationship is not as well established as the others, but the variation is probably 
regular enough to use for estimation of intermediate centerbody radii. 

Data for a specific boiler would be 


, boiler tube inside radius 

r , centerbody radius 
o 

p, pitch of the helical vane 

Other known data are p^, p^ , and C^. The latter constant may be taken as 0.4 
for a wide range of conditions. 


Fluid: Potassium at 1800 °F saturation temperature. 



Particular geometrical constants required are: 


& 

D 


D 

2r 


1 


4 


o 

Since the droplet starts its motion at r , and impacts on the tube wall at the radius 

o 

D/2, then r/r = 4 at impact. From Figure 38, t /t « 0.78 at impact. 

° r r i80° 

From Figure 37, 0 = 118° at impact. 

Consider droplets of the size such that 


D_ = 100 

D 

P 


100 






For the stated conditions, 


r 

o 



(108) 


= 0.15 

From Figure 39, t^ = 10 at impact. 

The axia\*normalized time to impact is, from Equation (107) , t 

£ 

From Figure 11, (for g negligible) z = 0.8. 

From the definition of z. 


=; 1.91. 


z z 



For the stated conditions, 

z , _ 

“ — . 833 z 


(109) 


The impact point will be at an axial distance down the tube corresponding to 


TT' .667 

An insert interruption only one tube diameter long would be adequate for this 
size of droplet. A graph showing the required spacing for other conditions is 
shown in Figure 40. It appears from these calculations that for the short pitch 
helical inserts generally proposed a one diameter long interruption would be 
adequate for even quite large droplets. 

A rather striking photograph of the deposition of droplets shed from the 
interrupted centerbody of a helical insert is shown in Figure 41. This photograph, 
which supports the analysis presented here, was taken during the air-water experi- 
ments reported in Reference 3, and was provided through the courtesy of Mr. G. 
Converse of General Electric Missile and Space Division. 
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Z/D , AXIAL DISTANCE TO IMPACT, IN DIAMETERS 



Figure 40. Axial distance to impact on wall as a function of droplet diameter 
helix pitch, and centerbody diameter for cases with negligible 
gravity field. 
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W 0 . = 0.00217 lb/sec Tube I.D. = 0.875 in. 

r^U 

W . = 0.050 lb/sec. x = 0.9584 

air 



Figure 41. Droplet deposition downstream from an interrupted centerbody. (Photo courtesy of 
General Electric Missile and Space Division, Cincinnati, Ohio.) 


11) The Motion of Very Small Droplets in Helical Flow 

The motion of large droplets in helical flow has been studied analyti- 
cally in Section B.2 through B.10. In these analyses the simplifying assumption 
was made that the velocity of the particle was small compared to that of the 
vapor flow. In Section B.10 it was shown that droplets large enough to satisfy 
that assumption are, for short p/D helical inserts and cases corresponding to 
actual boilers, deposited upon the tube wall in less than one tube diameter of axial length. 

The large slip analysis is not valid for all droplets. It is likely that 
a large portion of the droplets which might be entrained in the vapor flow are 
small enough so that starting from rest they attain the local vapor velocity in 
a distance which is short compared to the tube diameter. In such cases, the 
droplet velocities are almost equal to the vapor velocities throughout most of 
the flight time of the droplet. 


Such small droplets may be visualized as following the vapor flow 
but drifting toward the wall due to the radial force caused by helical flow. The 
radial motion through the vapor is opposed by drag and the droplet drifts without 
appreciable acceleration. The two radial forces may be equated and solved 
for the radial velocity as follows: 


F 


D 


H 


^ A 1 / dr \ _ v t 

D D2^v \dt / r 

where: 

F^ is the drag force exerted on the droplet by the vapor 

F is the radial force due to helical motion 
H 

is the droplet drag coefficient, taken as a constant 0.4 
Ap is the droplet frontal area 
Vp is the droplet volume 
is the vapor density 


( 110 ) 


1 


105 



is the liquid density 


r is the radius at which the droplet is situated at time t 
and v t is the rotational vapor velocity. 

For the idealized case of solid body rotation and uniform axial velocity, 


7 tD t n \ 

■ — V a(o r ) 


( 111 ) 


where 


and 


D is the tube inside diameter 

p is the helix pitch 

v is the axial mean velocity, 
a 


The simple differential Equation (110) maybe combined with (111), the result 
integi 
yield 


integrated with the boundary condition t = 0, r = r^, and rearranged to 



( 112 ) 


The droplet impacts on the wall when ^ 

r =_ 

The principal quantity of interest in this analysis is the axial distance of travel. 
Since it has been assumed that the slip is negligible and that the axial velocity 
is uniform, time is related to axial length by 


z = v & t (113) 

where z is distance measured along the tube axis. Equations (112) and (113) may 
be combined and rearranged to yield a relationship for the axial distance 
traveled by a droplet drifting from the tube centerline or centerbody radius to 
the wall. This result is 



( 114 ) 
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This relationship is plotted in Figure 42, using 


= 0.04 


which is characteristic of potassium at 1800 °F. 


The axial travel distance given by Equation (114) increases as the droplet size 
gets smaller . This is in contrast to the characteristic shown by the large 
slip analysis of Section BwlO. However, the distance to impact for both models 
is of the order of one tube diameter over a large range. 

Since the two analyses are based on two distinctly different idealiza- 
tions, it is important to establish ranges of applicability. It is reasonable to 
expect that the large slip model will prevail if the velocity of a droplet is not 
closely accomodated to tfyat of the vapor within one tube diameter. If the distance 
for velocity accomodation is short compared to one tube diameter, the small 
slip model will prevail. 

The relationship between axial location and velocity of a droplet introduced 
into a uniform flow is given by combination of Equations (51) and (52) of Section B.2 and is 


v 

P 



where v^ is the absolute droplet velocity. 

If adequate velocity accomodation is attained where 



0.98, then 


z 

D - 


45.1 


o D P 

L c JL 

4 D D p 
p * 


For potassium at 1800°F, 
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AXIAL DISTANCE TO I M PACT, IN DIAMETERS 






— = 1 when = 3.8 x 10^ 

D D 

P 

= 0.01 when = 3.8 x 10 5 

D D 

P 

-3 

One would expect the large slip analysis to hold for droplets larger than 10 D , 

-5 

and the small slip analysis to hold for droplets smaller than about 10 D . For 
droplet sizes in between these two values, the axial distance to impact should 
be bounded by the results of the two models. A composite graph of the two 
solutions is shown in Figure 43. 

In the small droplet small slip analysis only vapor rotation and vapor drag 
force are considered. However, these need not be the only mechanisms affect- 
ing drift of small particles. It might be supposed that the turbulent diffusion 
in the vapor flow, and the large secondary flows in helical systems (and other 
effects as well) would affect the transport of very small droplets toward the 
wall. 
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.AXIAL DISTANCE TO IMPACT ON WALL 






C. Liquid Distribution in Boilers 


1) Some Aspects of Two Phase Flow in Insert Transitions 

Some examples of the axial changes in boiler insert type which have 
been considered for use by designers are: fluted plug to helical vane, fluted 
plug to helical wire, and helical vane to helical wire. The vane and wire sec- 
tions may or may not have a center body associated with them, and the insert 
pitch may or may not be the same on both sides of the insert change. 

Most of these insert combinations would probably be used in regions 
in which the flow consists of two phases. The primary purpose of all of the 
inserts is to rotate the fluid, in order to separate the phases and to place the 
liquid component in contact with the heat transfer surface. For the purposes 
of this discussion it will be assumed that the upstream insert is very effec- 
tive in separating the phases, so that almost all of the liquid is in a continu- 
ous liquid layer on the wall. 

Consider the transition from a fluted plug to a helical vane system. 

If the pitches of both are the same, then the major difference between the 
two inserts is in the diameter of the center body, and the flow will be similar 
to a flow into a sudden enlargement. The vapor phase will decelerate rapidly, 
and there will be a recirculating stall at the base of the plug if the area change 
is large. The liquid layer on the wall will also tend to lose velocity and thick- 
en, although at a lesser rate than the vapor, and if the quality change with 
length is high enough no significant thickening of the liquid layer will occur. 

No observable effect of this transition upon the heat transfer performance 
would be expected. 

If the fluted plug and helical vane do not have the same pitch, the li- 
quid layer will be subjected to an acceleration which will thicken it on one side 


In 
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of the channel formed by the next insert and thin it on the other side. If the 
change in pitch is great enough, the liquid layer could separate from the vane 
entirely and leave a local dry spot on the heat transfer wall. This flow is iden- 
tical to the open channel flow around a sharp corner, except that the gravita- 
tional field is replaced by a centrifugal acceleration field which may be many 
times larger than gravity. 

From experience with open channel flows, one would expect that the 
major effects of the sharp corner turn would be obliterated within a couple of 
channel widths. If this transition were subjected to a high uniform heat flux 
a local hot spot might possibly develop at the dry spot; however, in a uniform 
wall temperature system, the dry spot would probably have no noticeable effect. 

One would expect the presence of a liquid layer on the tube wall to have 
a deleterious effect upon the operation of only one of the insert types listed a- 
bove, and that is the helical wire. Presuming that the liquid has by some means 
been deposited upon the wall, and presuming also that the liquid forms a layer 
which has appreciable thickness when compared to the diameter of the wire, 
the liquid layer will decrease the effective diameter of the wire. By the anal- 
ysis of Section A2, the force tending to rotate the fluid will be decreased by 
the liquid layer. If the liquid layer covers the wires, no force tending to ro- 
tate the fluid will exist at all. It is true that thick liquid layers will have a 
very high wave structure upon them, and the wires may be partially exposed 
even when the average liquid thickness is greater than the diameter of the wires, 
but based upon the analysis cited above, the force tending to rotate the vapor 
phase would be small. 

Consider the transition from a fluted plug or helical vane to a helical 
wire in the low quality region. The separation effectiveness of the plug and 
vane are both good so the liquid layer at the termination of these inserts may 
be thick. The effective transition the vapor phase undergoes in such a case 
could very well be from a helical vane (or plug) system to a very rough tube. 
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and one might expect the rotation to decay to virtually zero. 

The heat transfer in such a transition will not necessarily be poorer 
than that in a helical vane. As long as the wall is completely covered, liquid- 
like conductances will prevail. Since centrifugal accelerations are decreas- 
ing, more liquid may be entrained in the vapor, and surface waves will also 
be deeper. Both of these effects thin the liquid layer, and since a large por- 
tion of the vaporization in annular flow of liquid metals occurs at the interface 
between the liquid and vapor, a thin layer has less heat transfer resistance 
than a thick one. The irregular wave motion also increases the rate of tur- 
bulent diffusion of heat through the liquid layer. 

It appears that helical wires will be poor flow rotating inserts at low 
to intermediate vapor quality. The length of helical wire inserts in the low 
to intermediate quality region should be limited in length if downstream en- 
trainment of liquid in the vapor is expected to be a problem. 


2) Flow in a Curved Channel System 

An analysis of heat transfer and pressure drop in a curved channel sys- 
tem that represents the initial boiler tube insert (helical flow in an annulus) 
has been made. The results may be applied to the prediction of two phase pres- 
sure drop in the multichannel plug. 

The initial momentum transfer analysis is defined by separated, viscous 
liquid and viscous vapor flow (the liquid layer being held contiguous to the outer 
wall by centrifugal force). The classical shear stress-pressure gradient and 
viscous stress-strain equations for a curvilinear coordinate system are 


1 3p constant 

7* = — — + 

s 2 30 2 

r 


(117) 


and 



<U8) 


where t s shear stress 

O, angular position in the curved channel system 
r, radial position in the curved channel system 
v, curvilinear fluid velocity along a constant radius, r 
Equations (117) and (118) were solved utilizing four boundary conditions, namely, 
zero velocity at the inner and outer walls and the velocity and shear stress 
continuity at the gas- liquid interface. The resulting velocity profile solutions 
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across the liquid and vapor layers are shown below: 
for the outer, liquid layer, 


v = r 


for the inner, vapor layer, 
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A typical velocity profile together with the specific geometrical and fluid con- 
ditions is shown in Figure 44. This solution can be used to estimate the fric- 
tional pressure drop in a helical flow annulus (for a viscous-viscous flow). 
Addition of momentum and frictional pressure drop terms yield the total 
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pressure loss. The heat transfer for this system can be represented by con- 
duction through a superheated liquid layer or by the volume heat sink repre- 
sentation presented previously for annular flow in a tube (Reference 14 ). 

This same model is extended to the case in which the vapor phase is 
in turbulent flow in the following paragraphs. 

Integration of the shear stress-pressure and viscous shear stress- 
strain equations leads to the velocity equation for the viscous layer, namely, 

v = r ( ^ Inr - j ) + C 0 r (124) 

\2 U ]L B9 2 t i jL ) 2 

where and are constants determined .from the boundary conditions, r is 
radial position, n is the viscosity of the liquid, p is pressure and 6 is angular 
position. For the turbulent layer, the shear stress-pressure and shear stress- 
strain equations are again integrated, but in the latter equation this viscosity 
is replaced by an eddy diffusivity that is dependent on the mean velocity and 
radial position, namely, 

e = m v + B ( v ) (r A - r) (r - (125) 


where is the viscosity of the vapor, Bis a parameter that relates the tur- 
bulence level to the mean layer velocity and r^ is the radial distance to the li- 
quid vapor interface and r is the radius of the inner channel wall. The radius 
function in Equation 125 approximates the radial eddy diffusivity distribution. 
The previously indicated integration for the turbulent vapor layer velocity pro- 
file is 


(126) 
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The four constants, C^, C^> C^, and in Equations 124 and 126 are 

evaluated from the four boundary conditions, 

v (r 2 ) - 0 

v (r x ) = 0 

v^l) - v^") 

^ r i + ) = T J r P 

where r is the radius of the outer channel wall and r s is the shear stress. 

2 s 

The results are: 
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F 2 (r) 


F 3 (r) 



dr 
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dr 
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The integral in F (r) is of the form 
2 


F^(r) is of the form 


C dx 
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and the integral in 


(C 1 + °2 X + °3 x2> 


; they are expressible in long, known 


closed form algebraic equations available in standard integral tables. This solution 
can be used to estimate the frictional pressure loss in the initial boiler tube insert 
(helical flow in an annulus) . The solution also relates vapor quality to liquid layer 
thickness which can be used in a heat transfer analysis of the type described in 
Reference 14. 


3j A Criterion for the Onset of Entrainment Due to Vapor Shear at the Liquid- 
Vapor Interface in Annular Flow Boiling 

One of the important mechanisms causing the entrainment of liquid into the 

vapor phase from the liquid layer in annular flow boiling is shear entrainment. It 

appears that most of the shear entrainment is from the peaks of large "breaking" 

waves, which have a height above the average liquid layer thickness of the same order 

15 

as the layer thickness itself . There is no agreement as to the conditions upon 
which the onset of roll waves occurs. However, there has been a correlation published 
for the onset of spray entrainment from these waves which is purported to give good 

agreement with data for adiabatic two-fluid systems (probably air and water). This 

, . .12 

correlation is 
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where 


is the critical velocity for the onset of entrainment 

is gas (or vapor) phase viscosity 

is liquid surface tension 
is gas (or vapor) density 

is liquid density 

Although this correlation may give good results for the systems tested, 
and is claimed to represent systems with varied physical properties, it is easy 
to demonstrate that this correlation gives absurd results for some very impor- 
tant cases. 

Consider, for example, Freon 12 vapor and liquid at 155 psia. The 

critical velocity predicted by the correlation is 0.59 ft/sec. The dynamic head 

- 2 2 

for Freon 12 vapor at this velocity is about 2 x 10 lb/ ft , which is clearly 

far too low to entrain liquid from large waves. 

The corresponding predicted critical velocity for an air-water system 

at room temperature where one would expect the correlation to apply, is 78.5 ft/sec, 

2 

leading to a dynamic head of 10 lb/ ft , a reasonable value. 

It is possible to devise a simple model for the removal of liquid drop- 
lets from roll waves. Presume that the waves already exist on the liquid surface, 
and that the flow is upward. The removal of a droplet is probably due to the drag 
of the vapor upon the tip of an irregularity in the wave. A schematic drawing 


v 

c 



of this model is shown in Figure 45a. The force system acting upon the wave tip 
is indicated in Figure 46b. Clearly the restraining force on the wave tip decreases 
and the drag force increases rapidly after a "neck" appears, and a droplet is 
almost certain to be pulled off. A force balance may be written for the droplet 
when it is still attached to the wave tip. 

F_ = F + F (132) 

D g a 

The forces are approximately 


„ 7T 2 1 V 2 

F = C — D v 

D D 4 p 2 g v 

c 

where v is the local vapor velocity relative to the liquid 

f 4^t d o 3 

g g„ i 6 p 


F = 7r D a 

<y P 



Combination of Equations (1) and (2) and rearrangement yields 

1 

2V3" 

V v 3 

Evaluation of Equation (124) requires estimation of the droplet diameter and the 
drag coefficient for the wave tip. 

The droplet diameter is probably of the same order as the height of 
the wave tip above the mean liquid layer thickness. For a typical air-water 
system, the mean liquid layer must be at times in the order of 0.1 inch. The 
wave amplitude might be (as stated before) twice this value. So a droplet diameter 
of 0.1 inch would be a plausible value. 
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Figure 45. Liquid and vapor in upward flow. 
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The drag coefficient for a cylindrical object with a rounded downstream 
end is in the order of 0.25. 

Evaluating Equation (134) with these estimates and the other required 


constants 

cr = 5 x 10 3 lb /ft 
H 2° 1 

p =0.1 lb/ft 3 
v 

p„ = 62.4 lb/ft 3 
1 


The predicted value of critical air velocity is 


_ 1 _ 

2 

v = 1.15 1 670 + 4650 } ft/sec 
v 


= 84 ft/sec 


This is close to the value predicted by the experimental correlation. By use 

of the same droplet diameter and drag coefficient, a prediction may also be 

made for Freon 12 at 155 psia. The resulting critical velocity is 

_1 

v = 1.15127.5 + 13. 0} 2 
v 

= 6.8 ft/sec 

This is certainly a more reasonable number than before. The dynamic pres- 

2 

sure in the vapor based on this velocity is 2.3 lb f /ft . It is interesting that 
the gravity term which is small in the case of water with its high surface • 
tension, is the larger term in the case of Freon 12. 

The predicted entrainment criterion can be put into the same form 
as the experimental correlation. The result is 



6g 


1 _ 

~2 


up. C D 
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(135) 


122 



r 


The right hand side replaces the experimental constant in the correlation. 

Using the same physical constants as before, the right hand side has the value 

for air and water: _ _ -4 

2.8 x 10 

-*4. 

The experimentally determined value is 2.46 x 10 ' Such good agreement is 

probably fortuitous, since the model presented is simple and the phenomenon 
treated is quite complex. 

In the absence of a comprehensive theory of entrainment, the re- 
sults of this model have been used to make predictions for the boundaries of 
shear entrainment in a boiling system. These predictions are described in 
the next section. 

4) Estimated Limits for Shear Entrainment in Annular Flow Wetted- Wall Boiling Systems 

The mean thickness of the turbulent annular film in a forced convec- 
tion boiler may be predicted by means of mathematical models available in the 
technical literature (see for example Reference 14). 

The principal parameters upon which the liquid layer thickness de- 
pends are: vapor quality, vapor-to-liquid density ratio, wall friction factor, 
liquid-vapor interface friction factor, and tube diameter. 

The turbulent annular film may be expected to have a complex wave 
structure superimposed upon it, so that under some conditions large roll waves 
may occur and these large roll waves are apparently the source of liquid drop- 
lets entrained by shear. Since the "breaking" waves are of the order of twice 
the liquid layer thickness in height , it may be expected that the diameter of 
the sheared-off fragments is approximately the same as the mean liquid layer 
thickness. By combining the shear entrainment criterion presented in the pre- 
ceding section with the annular flow liquid film thickness predictions a map of 
conditions under which shear entrainment would be expected results. Several 
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such maps have been constructed. These are presented in Figures 46, 47, 48, 
and 49. 

As an example of the interpretation of these figures, consider the 
graph for Freon 12 flowing in a .45 inch tube. Figure 46. By hypothesis the 
film thickness must be at least equal to the droplet diameter for entrainment 
to occur. Therefore, a given size droplet will only be entrained into the vapor 
at qualities for which the film thickness is greater than that drop size. However, 
the vapor velocity must be high enough to entrain droplets of that size so the 
entrainment region lies above the line representing the entrainment criterion. 
Conditions for entrainment of 0.05 inch Freon droplets exist above the graph 
of the shear entrainment criterion, and below the 0.05 inch film thickness line, 
in the cross hatched area. The interpretation of the figures for the other fluids 
is the same except that they are complicated by lines representing the entrainment 
criterion for various sized droplets. By inspection of Figure 47, one would con- 
clude that shear entrainment will only occur for potassium in a 0.42 inch tube 
at low qualities and mass velocities far larger than those encountered in ex- 
perimental boilers. Figure 48 indicates that shear entrainment can occur more 
readily in a potassium flow in a 0.92 inch tube. Figure 49 which is constructed 
for wetting flow mercury indicates that shear entrainment is extremely unlikely 
under conditions of practical importance. 
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5) A Mechanism for Spray Generation in Nucleate Boiling 

There are several plausible mechanisms which may be proposed for 
the generation of spray by nucleate boiling. Careful observations of Freon 11 
in forced convection, water in pool boiling, and air bubbles rising in water 
have shown that one particular process occurs in a significant number of cases. 

In this process, bubbles rise to the surface and are caught in the surface layer 
where they form a blister. After remaining in the surface layer for a small 
length of time (up to ~ 10 seconds) they burst, and a small droplet is ejected 
from the surface. It seems reasonable to suppose that since the bubbles re- 
main in the surface layer for a finite length of time any upward momentum as- 
sociated with the rising bubble does not contribute significantly to the momentum 
of the ejected drop. 

The liquid layer covering the "blister 11 appears thinnest at the center, 
and it is reasonable to suppose that the surface ruptures at this point. One 
can therefore imagine a spherical cavity in the surface collapsing upon itself 
due to surface tension forces, and in the case of a gravitational field normal 
to the surface, due to energy in the ring-shaped wave formed by the broken 
bubble. When the bubble collapses upon itself, instantaneous pressures at the 
center must be high enough to eject a small amount of liquid from the surface. 

In the case of air bubbles in water, 1/8 inch diameter bubbles eject 

droplets to a height of two to three inches above the surface. These droplets 

-2 

appear to be in the order of 10 inches in diameter. 

Although no analytical description of the processes in the liquid layer 
following the fracture of the bubble seem to be presently available, the trajectory cf the 
liquid fragment ejected can be analyzed readily. 

Consider a spherical fragment leaving the surface vertically upward 

at the initial velocity v . The equation of motion of this particle is 

o 

dv 

" (F D + g e m) = m dt ( 136 > 

A first integration of this equation yields the following relationship between 
velocity and time: 


\i 
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where 


c i = — 
m v 


C 2 ge 


C 3 =V^C 2 tan_1 K V ° 


At the maximum altitude, v = o, so that 


t = C 

i max 3 


A second integration of this equation yields the following relationship between 
distance (from the surface) and time: 


• £ = C 4 - - log cosVC^ t d 


where 


‘d = S ' 1 


C 4 = 1 log cosve^ 2 C 3 


From these relationships, 


-i = log cos tan 1 
max 


°i 


It may be shown that in the case of the water droplet ejected vertically the fric- 
tional drag term is small compared to the gravity term. This is, of course, 
not true for a nucleation bubble in a vertically oriented film in forced convection. 
The constant C , is 

f d 

C 1 = — 
m v 


C p 
£ _D jv 
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P * 
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This will be recognized as the reciprocal of the characteristic length for drop- 
let acceleration used in several analyses in Section B of this report. 

For the water droplet observations described in this section, 

D - . 010 inch 
P 

Si =3 inches 
max 


so that 


C 


1 


0.036 

inches 


Jc C 
V 1 2 


3.74 

sec 



.965 x 10 


sec 

inch 


v * the initial velocity of the water droplet as it leaves the surface may be 
computed from these constants and Equation (139) . From this calculation: 


v =4.23 ft/sec 
o 

A droplet traveling upward three inches in a drag-free environment would have 
an initial velocity of 4.01 ft/sec. 
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6) A Discussion of the Mechanism by Which Liquid Collects on a Helical Flow 
Center Body 

A number of investigations, including recent ones at General Electric 

3 

performed under Contract NAS 3-9426, have shown that despite rather large 
radial accelerations because of rotational flow, a significant amount of liquid 
can collect on the center body of a helical flow system. Figure 50 clearly shows 
this center body flow. A qualitative description of the hydrodynamic causes for 
such a phenomenon follows. Included are suggestions for some simple modifications 
which can be made to the helical vane system to alleviate or completely eliminate 
such undesirable effects. A flow with curved streamlines has a pressure gradient 
normal to the streamlines . In a rotational flow system the flow streamlines are 
curved, and the pressure difference between the center body and the outside wall 
is given by r 

c 



A sketch of the geometry is given in Figure 51. This pressure difference 
would be uniform around the circumference of the tube if the velocity field were 
axially symmetric. There would also be no secondary flows. However, the 
vanes which maintain the flow rotation have boundary layers on either side and 
the value of the integral giving the pressure difference between the center body 
and outside wall tends toward zero as the vane is approached. There must 
therefore, be a pressure gradient component on the inside tube surface in the 
© direction, and strong secondary flows. 

An estimate of the pressure gradient on the tube wall in the 9 direction 
may be made. By consideration of Figure .52, 
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= 0.00685 lb/sec 


Tube I.D. = 0.875 in. 



W . = 0.0323 lb/sec 


C.B. = 0.25 in. 


x = 0.8250 


p/D =3.0 
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Figure 50. Two-phase flow in a tube containing a helical insert. Liauid layer on centerbody 
can be clearly seen. (Photo courtesy of General Electric Missile and Space 
Division, Cincinnati, Ohio.) 
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where n is the number of vanes on the helix. 

Clearly the pressure gradient gets larger as the number of vanes in- 
creases, and one would expect stronger secondary flows in multivane systems 
than in single vane systems. The secondary flows in tubes with helices have 
been investigated for single phase systems. As one might expect, the flow 
patterns are qualitatively as shown in Figure 53 „ 

Liquid droplets entrained in the vapor flow in a helical flow boiler are 
accelerated toward the outer wall and deposited there. The liquid flowing a- 
long the wall in a helical system has small velocity compared to the vapor and 
is propelled primarily by shear forces from the adjacent vapor flow. The 
boundary layer vapor flow, therefore, moves the fluid toward the vane and 
across the center body, and one observes a liquid rivulet in locations where 
transverse flows are low velocity; i. e. , on the center body at the stagnation 
region. In a boiler tube, the fluid which is not vaporized on the tube wall (the 
heat transfer surface) probably eventually migrates to the center body. Be- 
cause there is no heat generated on this surface, and there is no local force 
field strong enough to strip the fluid away from the center body, it simply 
fnigrates along the center body in the direction of the local shear forces. 

The flow toward the centerbody is difficult to show in photographs of a 
helical vane system. However, there is an insert system which has been tested 
by General Electric which demonstrates the flow toward the centerbody quite 
clearly. Figure 55 is a photograph of a tube containing helical wire flow rotator 
and a centerbody. These two objects are not in contact with each other, but a 
welded strut supports the centerbody at intervals. These struts provide a flow 
path for the liquid to cross over to the centerbody. The support strut near the 
left hand side of the photograph clearly has a stream of liquid flowing up it from 
the helical wire to the centerbody. 
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Figure 53. Secondary flow circulation pattern. 



Figure 54. Short single vane helix. 
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There are several obvious methods of alleviating this effect. One 
suggestion is to keep the number of vanes to a minimum (one) since in addition to 
causing the secondary flows, the vanes provide the path to the centerbody. A 
multichannel plug which has many paths to the centerbody available, would be 
expected to have a larger portion of the fluid on the centerbody than would a 
tube with a single vane. Observations of the SNAP-8 mercury boiler inlet plug, 
which has many small channels, have shown that the main portion of liquid 
leaving the plug is very effectively atomized. If this fluid were on the tube wall 
at the exit of the plug, which presumably is a region of fairly rapid vapor de- 
celleration (due to the increased flow area) one would not expect atomization 
to occur. On the other hand, a stream of relatively slowly moving liquid in- 
jected into the dividing line between a vapor jet and a strong recirculating eddy 
(which would be found 'at the base of the plug) would be atomized. One might, 
therefore, surmise that a large portion of the liquid leaving the SNAP-8 inlet 
plug may be on the center body rather than the tube wall. 

A method for further alleviating this effect would be to remove the center 
body entirely and terminate the vane short of the centerline, as shown in Fig- 
ure 54. The fluid would migrate to the inner edge of the vane, and probably 
form into droplets which would be stripped off and redeposited on the heat 
transfer surface. A large diameter helical wire would effectively approximate 
this geometry. Elimination of secondary flows would be accomplished by use 
of interrupted vanes, or propellers at intervals along the tube. A method of 
alleviating the collection of liquid on the center body while retaining the center 
body would be to put a skirt at intervals on the center body normal to the heli- 
cal vane. These skirts would force the rivulets to form droplets on the down- 
stream edge from which they would be stripped by the vapor flow. A sketch 
of the conceptual operation is shown in Figure 56. The possibility of using 
guide vane materials that do not make direct contact (wet) with the drops should 
also be considered; drop stripping by the vapor stream would probably be im- 
proved. 
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Figure 56. Drop-forming skirt for helix center body. 
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Motion of Wetted Liquid Films on Helical Vanes 


a) Secondary Vapor Flow Velocities 

It has been shown in the preceeding section that the longitudinal 
boundary layers on helical flow dividers provide the driving force for 
secondary flows. By consideration of the forces on a cross-section of the 
fluid, it is possible to estimate the velocities in the secondary flows. Such 
an estimate is made in the following paragraphs. 

Consider the idealized secondary flow depicted in Figure 57. The 
rotating fluid in the secondary flow has a center of motion about which it turns 
with a velocity profile approximating solid body rotation, but with a super- 
posed boundary layer against the solid surfaces. Idealizations necessary to 
perform the analysis are: 

1. The flow is incompressible , 

2. The frictional effects perpendicular to the tube axis and parallel to the 
tube axis are independent, 

3. Liquid layers which may be present on the walls have negligible velocity, 

4. Boundary layers in the tube may be approximated by one-seventh power 
velocity profiles, 

5. The secondary flows have no effect upon the idealized solid body 
rotation (with boundary layers) of the mean flow about the tube axis, 
which is caused by the helical divider. 
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Figure 57. Idealized flow patterns 



Figure 58. Axial velocity c 
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along line a-b. 


where 


and 


r^ is the tube radius 

is the velocity component about the tube axis 

r is radius from the centerline 

p is vapor density 


with the solid body rotation idealization, 


v t (r , y) 
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where v^(y) is the local axial velocity. 


(143) 


If one visualizes a conventional one-seventh power boundary layer to 
exist upon the helical divider, as shown in Figure 58, then 
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where y is distance from the divider 

D is the passage equivalent diameter for turbulent helical flow 


v is the maximum velocity (axial) 
max 


By combination of equations (142) through (144) and integration, the 

* 

radial pressure difference may be shown to be (neglecting the effect of 
boundary layer on the tube wall) 
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The integrations are not truly radial, but are taken parallel to the helical divider, 
which lies along a tube diameter. Geometrical inaccuracies such as this are 
probably negligible compared to other simplifications. 
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The mean axial velocity component is 


where 


and 


w 

v 



w^ is the vapor phase flow rate 
is vapor density 

A is the flow cross-section area. 
T 


The mean and maximum velocities are simply related 


11 


v = .817 v 
a max 


( 146 ) 
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One may then write 
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This equation describes the pressure along the tube wall extending away from 
point c , as it is affected by the axial boundary layer along the divider between 
points a and c. 


The pressure integral along the line a - b, where there is no boundary layer 
in the axial flow, is 


^ P a-b _ / 1 \ 2 / ttD \ 2 
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The usual equation used by most authors which neglects tube boundary layers is 
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(150) 


Although this latter equation gives excellent results, the form 
Equation (149) containing the correction for maximum velocity with a boundary 
layer is used here to be consistent with other equations in the analysis. 
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The boundary layer on the vane causes a pressure defect along the tube wall, 
and this pressure defect may be integrated along the wall away from the vane 
to obtain a net radial inward force being exerted on the fluid. The magnitude 
of this force driving the secondary flows per unit of tube axial length, is 
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This force, with the appropriate moment arm around the center of secondary 
flow rotation, provides a torque on the fluid which is restrained by the 
secondary flow shear stresses along the tube wall and helical vane (see Figure 59) 
and by turning losses due to eddying in the corners. The restraining forces per 
unit of tube length are in the order of : 
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where v is the velocity of the secondary flow, measured normal to the tube axis 
s 

just outside of the secondary flow boundary layer. 
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NO SHEAR ALONG 
THIS SEGMENT 


Figure 59. Transverse velocity profiles in secondary 


flow cell. 


The turning losses are taken to have a loss coefficient of one, there are 

effectively three corners, and the effective width of the turning stream is 

estimated at D 
e_ 

4 

For equilibrium secondary flows, these forces may be equated (since the 
moment arms are of equal length )and solved for the secondary flow velocity. 

These operations yield 
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characteristic of tight twisted dividers , a typical value is 
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for a typical £ value of 0.02 the denominator bracket terms are 
Jo. 06 + 1.5j 


Equation (154) reduces to 
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This simple relationship is based upon what is considered to be a reasonable 
flow model, so that although the coefficient is only an estimate, the functional 
form is probably correct to the first order. 

16 

Seymour has made extensive measurements of secondary flows in systems 

containing twisted tapes of quite long pitch. For a tape of ~ = 0.105, 

P 

his measurement yielded v = 25 ft/sec, v = 208 ft/sec, so that 




experimental 


0.12 


From Equation (155), the predicted value would be 


analytical 

The agreement is reasonable. 

The predicted pressure rise (from centerline to wall), computed from 

Equation (150) ,does not include the effects of boundary layers. Boundary 

layers would be expected to increase the measured pressure rise by increasing 

the maximum velocity and by decreasing the tube radius at which the maximum 

velocity occurs. Why this is so may be seen from Equation (142). It is 

worthwhile to assess the magnitude of this effect. For this calculation, 

suppose the axial component to be modified by a one-seventh power velocity 

profile on all surfaces. A mass flow integration over the half tube has been 

performed numerically to determine v . The result of this calculation is: 

max 
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1.33. 


v 

max 


v 

a 


This compares with the full tube value 


v . 

max _ _1 

v ~ .817 

a 


= 1.22 


Using the hypothetical one-seventh power profile with a maximum velocity of 

1.33 v , the pressure rise along a radius was computed numerically, by use 
a 

of Equations (142) and (143). The result was: 



boundary layers 
no boundary layers 


1.13 


It appears that differences between the ideal relationship Equation ( 150 ) and 
actual measurements will not be very large. 


b) Force Systems Acting on the Liquid Film 

Consider a wetted liquid film on a helical vane (or twisted tape) in a tube with 

a vapor flow. The film is driven in motion by shear from the vapor flow. 

This shear has two components, one at constant radius from the tube centerline 

and directed downstream, and one directed inward along the tube radius. The 
second component is due to the secondary flows in the tube. The liquid film 
moves along the vane due to the axial shear, and as it moves it rotates about 
the tube axis, so that there is a centrifugal force on the liquid tending to 
accelerate it toward the tube wall. Both radial and axial motion are opposed 
by friction between the vane and the liquid film. A schematic representation 
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of this force system is shown in Figure 60. The purpose of this analysis is 
to determine under what conditions the liquid film will flow toward the center- 
body. The flow and force systems in the axial and radial directions are 
considered to be independent of each other. 


A boundary layer analysis of the flow in the direction of constant radius 
yields (see, for example, reference 17 ) for a laminar film. 
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(156) 


where 

v is the velocity in the film in the direction of constant radius 
x 

is the vapor shear force in the x direction 

y is the distance from the vane surface (within the film) 
p ^ is liquid viscosity 


The radial acceleration to which the film is subjected is 
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a = 
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The radial flow equation, which includes this body force, is 
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Figure 60. Forces acting on an element of the liquid layer 



where 


is velocity in the radial direction 
p is liquid density 

The appropriate boundary conditions are 

y = 6 : 


and 
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y = 0 

v = 0 
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where 

<5 is the film thickness 

F is the radial shear force due to secondary flows. 


By integrating twice, and applying the boundary conditions, one obtains 
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The quantity of interest is the net flow of fluid in the radial direction. This 
flow, per unit of vane length, is 
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The radial flow outward goes negative when the second term is greater 
9 

than ~ , and the net flow is toward the centerbody. 


The vapor shear stress due to secondary flows is (from part a) of this section) 
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It is clear from inspection that there is always a value of 6 , the 
film thickness, small enough so that Equation (165) will be greater than 
9/60, regardless of the values of the other parameters. 

The net-flow of a "thick" film is outward, the net flow of a "thin" 
film is inward. One may conclude that when the liquid wets the wall, there 
will always be a net flow of liquid toward the centerbody in a two-phase 
helical flow system. 

8 ) Criteria for the Existence of the Liquid Phase on Boiler Tube Centerbodies 

The criteria for the existence of liquid phase on centerbodies in boiler 
tubes can best be outlined by grouping the processes into hypothetical steps. 

First, it is postulated that a liquid layer exists on the centerbody as 
well as on the boiler tube wall. In this case vapor flows between the two 
liquid layers. Linear and curved channel three-layer flow analyses similar 
to the two-layer solutions previously discussed can be made. From resulting 
liquid film thicknesses on the centerbody, tests on liquid film rupture can be 
made and the fraction of centerbody area covered by rivulets determined 
(following the procedure used for the boiler tube wall) . In the case of curvilinear 
flow, however, the previous analysis for the rivulet profile would have to be 
generalized. 

Specifically, a centrifugal force as a result of rivulet flow around the 
centerbody would be added to the force equation. As the centrifugal force term 
increases in strength, the rivulet profile will become thicker and eventually 
the centrifugal force will exceed the surface tension force and the rivulet will 
leave the surface or rupture. This simple criterion suggests practical ways of 
reducing the liquid flow on centerbodies. If the contact angle for a given liquid 
metal and surface combination were to be increased by choosing an appropriate 
material for the centerbody, rivulets would be removed from the surface or their 
density decreased at reduced vapor flow rates or qualities. 
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D. Temperature Distributions in Boiler Tube Walls 

It is of technical interest to be able to describe quantitatively the 
peripheral variation in wall temperature in a boiler tube that contains 
discontinuous liquid layers. These layers may be in direct contact with the 
wall or be separated from the wall by a vapor film. The difference in heat 
transfer between the liquid cooled zone and the vapor cooled zone contiguous 
to it causes a temperature asymmetry. An analytical solution has been derived 
for a thin-walled system; the details and results are given below. Analog 
solutions for temperature distributions in thick-walled boiler tubes subjected to 
several liquid distributions are presented later in this section. 

1. Temperature Distributions in a Thin-Walled Tube with Internal 

Conductance Discontinuities 

The idealized boiler tube system is shown in Figure 61 and the postulates 
that define the system follow: 

a) The boiler tube wall is thin; thus, the temperature drop across 
the wall is small and negligible compared to peripheral temperature drops. 

b) The ends of the one-dimensional wall element shown in Figure 61 
are insulated. 

c) A step function conductance profile exists on the boiling side, 
simulating the rivulet and vapor zones. 

d) A uniform conductance distribution exists on the hot fluid side. 

e) Steady state exists. 

A heat rate balance on a differential wall element (which accounts for 
convection and conduction into and out of the element) yields, 
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The solutions of Equations (166) and (167) are, respectively, 
for 0 < x < x 
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C r C . C 3 , c , equation constants (see derivation) 

2 

convective conductance on hot fluid side, Btu/hr ft °F 

2 

h y > vapor convective conductance, Btu/hr ft °F 

2 

h £ , rivulet vaporization conductance, Btu/hr ft °F 
x 

k .wall thermal conductivity, Btu/hr ft °F 
w 

n , n , a , a , system parameters (see derivation) 
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T *boiler wall temperature at position x, °F 
w 

x , peripheral distance along boiler tube (origin at rivulet center) , ft 

x , half width of rivulet, ft 
r 

x , peripheral distance from rivulet center to point where the peripheral tem- 
perature gradient is zero, ft 
6, boiler tube wall thickness, ft 

T,, mixed mean temperature of the hot fluid, °F 
and “ 

T , saturated vapor temperature, °F 
s 

dT 

From the boundary condition, ^ = 0 at x = 0, it can be shown that 



dT 

From the boundary condition, — ■ - = 0 at x = x , it can be shown that 
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equations in two unknowns (C and C ) result, 
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The solution of Equations (172) and (173) yields 
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Thus the solution to the boundary value problem described above is given by 
Equations (168), (169), (170), (171), (174), (175), and (176). 
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The solution has been evaluated for three specific examples. 

The results of the three sets of calculations are shown in Figures 62, 63, and 64 
to demonstrate how the shapes of the wall temperature profiles can vary as 
the system parameters are changed. The thermal conductivities of the boiler 
tube wall range from 0 to « and fractions of the wall covered by a rivulet of 0. 1 

and 0.5 are considered. The vapor conductance in the boiler tube, h , varies 

2 v 
from 10 to 100 Btu/hr ft °F. The hot fluid T, is constant at 1300°F and the 

h 

saturated vapor temperature in the boiler tube, T , is constant at 1000°F. The 

s 

hot fluid conductance, h^, and the rivulet conductance in the boiler tube, h^ , 
are both set equal to 10,000 Btu/hr ft^ °F. The boiler tube wall thickness, <5 , 
is equal to 0.004 ft. The specific variations considered in Figures 62,; 63, 
and 64, are given in the legends. 

2. Temperature Distributions Due to Axial Discontinuities in the 
Internal Liquid Film in a Thick-walled Boiler 

Boilers which are wet by the vaporizing liquid in the low quality region 
all undergo a transition at some quality in which the wall ceases to be completely 
wetted by liquid. If this point occurs at a quality of less than one, the remain- 
ing liquid may be in rivulets on the wall, (which is locally wet) rivulets separated 
from the wall by a vapor film, droplets or masses of liquid throughout the bulk 
of the vapor, or combinations of these distributions. 

Consider the point at which the liquid film ends, so that the wall cooling 
is discontinuous on the inner surface. If the wall conductivity in the direction 
normal to the discontinuity is significant then heat will flow along the wall toward 
the high conductance (low temperature) region, increasing the local temperature 
difference in the well cooled portion near the discontinuity. The magnitude of the 
inner wall heat flux and temperature variations are of interest. 

The solution for the temperature and heat flux distribution in a thermal 
conductor subjected to the appropriate boundary conditions is not easily obtained 
by analytical means, and the required information was obtained by a two- 
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TEMPERATURE, °F 



Figure 62. Temperature profiles foi a thin wall boiler tube with a rivulet, 
(X p = 0.01', X e = 0.02', h y = 100 Btu/hr ft 2 °F.) 
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dimensional electric analog, using Teledeltos (Western Electric Company) con- 
ductive paper. The circuit used is shown in Figure 65. A photograph of the ap- 
paratus is shown in Figure 66. A tube cross section model for study of tempera- 
ture under rivulets, which is the case discussed in the next section, is shown 
in the photograph. 

All of the model work was done using the uniform outside wall heat flux 
boundary condition because this condition is frequently used in 'critical heat 
flux* boiling experiments and also because it eliminates the outside wall con- 
ductance as a parameter. 

The two parameters used are. 





which is the ratio of the 'wet wall* conductance to the reciprocal of boiler wall 
resistance and 



h 

v 


which is the ratio of the 'wet wall' conductance to the 'dry wall' conductance, 
h is the heat transfer conductance, k w is the wall thermal conductivity and <5 
is the wall thickness. 

Two typical wall temperature distribution plots are shown in Figures 67 and 
68. If is clear from these plots that although the outside wall heat flux is uniform, 
the inside wall heat flux is highly non-uniform for many wall thicknesses in 
either direction from the discontinuity. At large distance from the discontinuity 
the isotherms lie parallel to the wall edges. In the portions of the models 
shown, the isotherms lie almost perpendicular to the wall edges, showing that 
the major portion of the heat is flowing along the wall. In Figure 69 the variation 
in inside wall temperature at the conductance discontinuity is shown as a function 
of the two parameters. 

These results may be related to an actual boiler tube under uniform heat 
flux conditions as follows. The wall temperature drop at x = ± °° is 
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Figure 65. Circuit for flux plotting apparatus. 
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Figure 67 . Temperature distribution in wall near conductance discontinuity 
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Figure 69. Normalized temperatui 
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at conductance discontinuity. 




The temperature of the inside wall at - 00 (the 'wet wall' end) is given by, 


(T 


inside wall at - 00 


) 



The temperature T is given by 

T, -4T, (q ^ A) + l l' ,A + T 

> \ 

where T g is the saturation temperature. 
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3 • Temperature Distributions in a Thick-Walled Tube with Internal 
Rivulet Flows 

a) Number of Rivulets Variable 

A plausible internal phase distribution can be visualized in which the 
axial variation in inside conductance is negligible, but in which the azimuthal 
variation is large, due to wetted liquid rivulets running along the inner wall of 
the tube. In general, the significant parameters in this hypothetical flow dis- 
tribution are: boiler wall thickness 6 , wall thermal conductivity k^, wetted 
wall conductance h^ which may be a function of location in the rivulet, 

’dry wall' conductance h^ which may be a function also of azimuthal posi- 
tion, inside to outside tube radius ratio, ratio of rivulet width to inside cir- 
cumference, number of rivulets, and probably others as well. The general 
case contains too many parameters to present concisely, and a fairly specific 
case has been investigated. General Electric Company SPPS performed ex- 
periments with a 3/8 schedule 80 (IPS) columbium - 1% zirconium boiler tube under 
NASA contract NAS 3-2528, and this boiler tube was chosen to model with the 
flux plots. Radius ratio, wall thickness, wall thermal conductivity are all fixed. 

In addition, fixed representative values of 'wetted wall' conductance and 'dry 
wall' conductance were chosen. The inside conductances were idealized as 
uniform within their respective areas. A single rivulet was taken to be 1/10 
of the inside circumference. The only variable then considered was the number 
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of rivulets on the inside wall- The fixed parameters are 

tube radius ratio = 0.63 
6 = 0. 126 inch 


Btu 

k w = 31,6 ft 2 hr CF/ft) 


= 12,000 


Btu 

ft 2 hr °F 


h 

v 


= 600 


Btu 

ft 2 hr °F 


Based on these values, the parameters used in the flux plot work were taken as 



k/6 


4 


= 20 


v 

Flux plots were made for cases of 1, 2, 3, and 5 idealized rivulets. The case in 
which all the wall is wetted except for one dry strip 1/10 of the inside circum- 
ference wide was also plotted. 

The 5 flux plots resulting from this work are given in Figures 70, 72, 74, 
76, and 78 of this report. The outside and inside wall temperature variations 
for these same cases are given in Figures 71, 73, 75, 77, and 79. In all cases 
the temperature distributions have been normalized upon the temperature drop 
across the wall which would have existed had the inside wall cooling been uniform. 

Inside wall to saturated fluid temperature differences for the center and 
edge of the rivulet are given in Figure 80. These temperatures have been nor- 
malized by dividing by the temperature difference which would have existed had 
the entire inside of the tube been wetted. Outside wall maximum and minimum 
temperatures are given in Figure 81. 
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Although this model is rather specific several general conclusions can 
be drawn from the results . One conclusion is that when wetted rivulet flow 
exists in a boiler tube, very large scatter in the conductance data could occur 
if the conductances are computed from a single outside thermocouple. A sec- 
ond conclusion is that onset of film boiling of the rivulets could occur at 
smaller apparent inside wall to fluid temperature differences than one might 
expect. Under conditions of non-uniform cooling the boiler tube may be sub- 
jected to considerable thermal stresses and transverse warping might also 
occur. 
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Figure 74. Boiler surface temperatures. 
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Figure 75. Boiler surface temperatures. 
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Figure 79. Boiler surface temperatures. 
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Figure 80. 


Local inside wall to saturated fluid temperature 
differences, normalized on uniform "wet wall" 
cooling temperature difference. 
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Figure 81. Outside wall temperature maxima and minima 



b) Width of Rivulets Variable 

In the preceeding section surface and interior wall temperature 
distributions were presented for cases in which the maximum and minimum inside 
conductances and hypothetical rivulet width were held constant. The number of 
rivulets was the only variable. Since general results cannot be obtained in a problem 
with as many variables as are present in rivulet flow heat transfer, another specific 
case is presented here for comparison with the preceeding one in order to test the 
generality of the conclusions. In the present case, the geometry and conductances 
are identical to those in Section C . 3 . a) , but the number of rivulets is held at two 
and the width of the rivulets is allowed to vary. 

The normalized wall temperature distributions on the inside and 
outside surfaces are shown in Figure 82. The variation of maximum and minimum 
wall temperatures with amount of inside area Wetted are shown in Figures 83 and 84. 

In Figure 81, the number of rivulets corresponds to the percentage 
area covered by liquid, so that two rivulets, for example, correspond to 20% wetted 
area or 80% dry. Then by comparing Figures 81 and 83, it may be seen that the 
variations in outside wall maxima and minima are similar for the two cases, and 
it may be possible to generalize the conclusions to arbitrary liquid rivulet distributions. 
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Figure 82. Outside and insid 





Figure 83. Outside wall maximum and minimum temperatures 








4) Apparent Heat Transfer Conductances in a Boiler Tube Containing Rivulet Flows 

Heat transfer conductances measured for boiling flows are usually based on 

data obtained at the outside wall of the boiler tube, so that both the heat flux 

and wall temperature calculations for the inside wall require a correction for 

changes across the wall. For lack of other information, this correction is 

based on the assumption that the heat flow is purely radial. It has been shown 

in the preceeding sections that in the presence of rivulet flow, or near the termination of 

annular flow, large variations in heat flux distribution can occur on the inside 
surface even if the flux on the outside is precisely uniform. The conductances 
inferred are therefore incorrect and it is useful to assess the magnitude of this 
error. 

Consider a boiler tube containing a rivulet flow and subjected to a 
uniform outside heat flux. The rivulets are considered to be axially uniform so 
that all variations are two dimensional. The conductance between the rivulet 
and the wall is h^ , the conductance between the vapor and the wall is h y , and 
the portions of the inside tube area covered by the rivulets and vapor are A ^ and 
A^ respectively. The apparent conductance on the inside wall, based on outside 
wall data, is 

(q/A) o (A q /A.) ( 180 ) 

h a = T : _ T - AT 

o s w 
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where (q/ A) is the uniform outside heat flux 

o 

A is outside wall area 
o 

A^ is inside wall area 

T is outside wall temperature (this varies with angular position and is 
assumed to be measured from a single, arbitrarily placed 
thermocouple) 

T g is saturation temperature 
AT^ is wall temperature drop 

For usual boiler tube geometries, AT^ is probably adequately determined by 


AT = 
w 



( 181 ) 


where k is the wall thermal conductivity 
6 is the wall thickness 


The actual mean conductance on the inner wall is 
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or, because of the known relationships between 
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For cases in which the rivulet conductance is much larger than the vapor con- 


ductance, this relationship reduces to 
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The maximum and minimum possible apparent conductances have been computed 
from the wall temperature distributions of Section c.3.b), using Equations (180) 
and (181). These have been normalized using Equation (183), and are presented in 
Figure 85. It may be seen that the possible scatter band in the apparent conduct- 
ance with this thick boiler tube is only about 40% wide and a randomly placed 
thermocouple would be expected to yield a conductance below the true mean value. 

The maximum and minimum conductances based on the data of 
Section C.3.a) for constant width rivulets of varying number, are shown in Figure 86. 
A continuous curve is not valid for this case, but it may be seen that the maximum 
expected error is about the same and again the randomly placed thermocouple 
would usually yield a low conductance. 

One would expect that scatter in mean conductance would be increased 
for thinner walled boilers in spite of the smaller heat conduction around the wall. 

For a vanishingly thin wall, where the heat flow is purely radial, the scatter in 
mean conductance inferred from a randomly placed thermocouple would be as 
large as the maximum and minimum conductances themselves. 
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Figure 85. Error in rivulet flow boiling conductance inferred 
from outside wall data, columbium boiler. 
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Figure 86. Error in rivulet flow boiling conductance inferred 
from outside wall data, columbium boiler. 
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E. Mechanisms of Heat Transfer Regime Change 

Analyses and experimental data are available for prediction of heat transfer 
in many of the flow regimes known to exist in forced convection boilers. It is not 
always clear however, which of the many possible flow patterns will occur at a 
given point in a boiler. The analyses in this section have been worked out in order 
to provide some information about which flow regime might be expected under a 
given set of circumstances. 

1) Hydrodynamics of Linear Wetting Flow Boiling "Burnout" 

The term "burnout" , as it is usually used in reference to forced con- 
vection boiling of wetting fluids in vertical linear flow systems, is assumed 
by most workers to coincide with the end of the annular liquid film. Obser- 
vations made with Freon 11 in a glass test sfection at Geoscience Limited 
(Reference 14) support this interpretation of "burnout". These observations, 
taken under uniform heat flux conditions, demonstrate a very large and dis- 
continuous decrease in heat transfer conductance at an axial location which 
coincides with the termination of the annular film. The "burnout heat flux" 
data obtained from these experiments coincide with the experimental critical 
heat flux data reported by Kirby (Reference 12) for Freon 12. These latter' 
data are reproduced in Figure 87. One may conclude that in Kirby’s data, 
which were taken in an opaque test section, the "burnout quality" coincides 
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with the downstream termination of the liquid film. 

As has been pointed out by Kirby, the rate of change of the liquid flow 
in the film is determined by a mass balance between the liquid evaporated, the 
entrained liquid deposited from the vapor flow, and the film liquid entrained 
into the vapor flow. The liquid entrained from the film is usually supposed 
to be removed from the film by two mechanisms; namely stripping of liquid 
from wave peaks, and ejection of liquid by breaking of bubbles generated in 
the film due to boiling. If such a picture of the annular film is actually true, 
then the mechanisms controlling the flow rate in the annular film are largely 
hydrodynamic. One might expect that the liquid film flow rate at any point 
in the boiler is determined by the processes occurring upstream; the rate 
of deposition, the rate of entrainment, the rate of change of vapor velocity 
all contribute to the liquid film flow. One might also expect that there are 
many different combinations of processes which would give virtually the same 
liquid flow rate at a particular point. 

The annular film is very thin near its end so that bubble nucleation 
and wave motion are probably both suppressed. If such is the case, then the 
principal mass fluxes in and out of the annular film near its termination are 
due to evaporation and deposition. When the end of the annular film occurs at 
a quality less than one, the remainder of the liquid must certainly be entrained 
in the vapor flow. Due to diffusion, the fine droplets in the vapor eventually 
settle out on the walls and evaporate. The longer is the boiler, the smaller 
the entrained flow at the end will be. This deposition of entrained liquid from 
the vapor flow onto the walls would not be expected to be greatly affected by 
the presence of a thin liquid film on the wall, so that the deposition process 
may be considered without regard to the presence of a liquid film. 

In a uniform heat flux boiler, length, quality, and heat flux bear a sim- 
ple relationship to each other, namely 
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( 184 ) 



It is clear that for a fixed flow rate, a low heat flux boiler is a long boiler. 

If entrainment rates are not too high, a long boiler Las proportionately less 
entrained liquid at a ^ . -a quality, and one would expect that the end of the 
annular fi 1 ' .,^aid occur at higher quality for a low flux (long) boiler than for 
a high flux (short) boiler. 

Some conclusions may be drawn concerning the nature of "burnout 
heAt flux" experimental determinations in linear flow. One is that variations 
upon the conventional straight, vertical) continuously-heated tube geometry 
may be expected to change the deposition and entrainment rates. These vari- 
ations will therefore change the resulting "critical" heat flux. Thus a curved, 
or unheated section in the two phase region would be expected to decrease the 
amount of entrained fluid, thicken the film, modify the entrainment rate, and 
change the quality at which the end of the annular film occurs. Another con- 
clusion is that regardless of tube diameter, liquid flow rate, or other condi- 
tions, the end of the annular film will occur atX = 1 for vanishingly small 
heat flux. 

Some of the observations made above may seem rather obvious. jIow- 

ever, they are appropriate to use as a guide for the interpretation of "buri out" 

data. For example, the plot of Figure 87, showing Freon 12 burnout data, 

indicates a linear relationship between (q/A) andx over a wide range of 

B B 

test conditions. From the physical arguments just given, it may be concluded 
that all of the data lines must pass through the point (q/A) = 0, X = 1. 

Since this is the case, only one of the lines in Figure 87 can be straight at sub- 
stantially lower flux than the data cover. All of the remainder must curve at 
low heat flux. Since this is the case, extrapolation of these data to lower fluxes 
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than the range of the data is completely unwarranted. The data of Figure 87, 
plotted on semi-logarithmic coordinates, are shown in Figure 88. It is clearer 
in this plot that the data do tend toward zero heat flux at X = 1 . The solid lines 
indicate the actual range of the data. The dotted lines indicate where one might 
expect data taken over a wider range of heat fluxes to lie. These same data 
have been replotted as functions of length by use of Equation (184) , and are pre- 
sented in Figure 89. In this figure it may be seen that the boiler length re- 
quired to maintain the annular film at high quality increases rapidly for high 
mass flow rates. 

One of the unexpected aspects of these burnout data is the amount of 
liquid entrained in the vapor flow at high mass flow rates. As stated earlier, 
all of the remaining liquid is entrained at the end of the film. Thus, the en- 
trained liquid plot of Figure 90 has been constructed from the data of Figure 89 
by use of the relationship 


W 

e 


W. 



(185) 


where W is the entrained liquid flow rate and W, is the total flow. From Fig- 
e t 

ure 90 it can be seen that for high flows and fluxes as much as 90% of the flow 
can be entrained. When studying these data, it is well to keep in mind that these 
graphs of entrained liquid do not represent variation of entrainment with length 
in a single boiler at a given condition, but represent instead the entrainment at 
the end of the annular film in tubes of fixed diameter and flow rate but of many 
different lengths (heat fluxes). 
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Figure 88 . 


Freon 12 linear flow "burnout" data. (From Reference 12) 














































2. A Criterion for the Onset of Nucleate Boiling in a Forced Convection System 
With an Annular Liquid Film 

The phenomenon of nucleate boiling in conventional pool boiling systems 
is usually attributed to minute cavities with residual gas trapped inside. These 
minute cavities form the visible nucleation sites which are the source of vapor 
bubbles in the pool. Water usually boils in metal containers with only a degree 
or two of superheat at the wall. However, it has the ability to superheat in the 
liquid phase to almost 150 °F above the saturation temperature if carefully de- 
gassed and heated in a system which has been evacuated to remove the trapped 
gas from the cavities in the heated surface. In contrast to metal, glass usually 
has few cavities in the surface which can act as nucleation sites. Boiling chips, 
which contain many microscopic cavities, are frequently added to glass boiling 
systems to reduce liquid superheating with its accompanying violent change of 
phase. 

Consider the boiling process in a forced convection system containing 
an annular film. If nucleate boiling does not occur then heat which is added to 
the annular film at the wall must be conducted across the liquid layer so that 
vaporization occurs on the inner surface at the interface between the liquid and 
the vapor. Since this interface must be at or above the local saturation tempera- 
ture, the liquid in contact with the heated wall must be superheated. One would 
expect that in the absence of nucleation sites, the onset of nucleate boiling would 
be dependent upon the maximum superheat temperature which the fluid can sus- 
tain without changing phase. In the forced convection experiments of Reference 18, 
Freon 22 was vaporized in three tubes of identical geometry and under very nearly 
identical conditions, except that two of the test sections were metal and the third 
was glass. The glass test section was heated by electrical conduction through a 
transparent metal film on the inside and was used for flow visualization. Some 
of the data from these experiments are reproduced in Figures 91 and 92. Clearly 
in these data the glass has a much lower boiling heat transfer conductance than 
the two metal tubes. Since this difference must be related to the surface material, 
it is a plausible conclusion that the difference is in the availability of nucleation 
sites in the metal and their absence in the glass test section. Since the glass test 
section by hypothesis does not have as many nucleation sites as the metal tube, 
more superheat of the liquid at the wall would be required to bring on vapor genera- 
tion in the liquid bulk. 
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Figure 91. (Reproduced from Figure 8, Referencel 
-22 at a reduced pressure of 0.12. R; 
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) Heat transfer data for refrigerant 
ige of vapor qualities from 0.01 to 



LOCAL HEAT FLUX ( BTU / HR FT^) 


VERTICAL STAINLESS STEEL T SAT =82"F 
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Figure 92. (Reproduced from Figure 9. Reference !8) Heat transfer data for refrigerent 
-22 at a reduced pressure of 0.22. Range of vapor qualities from 0.01 to 
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they are charged with the working fluid. This procedure is in contrast to that 
used in conventional fluid systems. The evacuation is required in order to pre- 
vent chemical reactions or oxidation of the working fluid. Since the removal of 
inert or noncondensible gases from the system by evacuation may destroy the 
effectiveness of small surface cavities as nucleation sites, one would expect 
that superheating of the liquid metal would occur readily in liquid metal systems. 
Many workers have reported and described this phenomenon. 

In the complete absence of nucleation sites it is possible to estimate 
the amount of superheat required to bring on the formation of vapor cavities in 
the liquid. Predictions have been made for water, potassium, and sodium in 
Reference 19. These results are reproduced in Figure 93. In the absence of 
nucleate boiling, (which acts as a volume heat sink within the liquid film), all 
of the heat added to the annular film at the wall must be conducted or convected 
across the liquid layer to vaporize fluid at the inner surface. Since the liquid 
metals have very high thermal conductivity, convection may be neglected and 
the temperature drop across the liquid layer computed on the basis of conduc- 
tion alone. With this simple model, graphs for the minimum heat flux required 
to bring on nucleate boiling have been computed for potassium and Freon 12 
and are shown in Figures 94 and 95 assuming that the interface between liquid 
layer and vapor core is at the saturation temperature. The curve for Freon is, of 
course, hypothetical since with its low thermal conductivity turbulent convection 
in the liquid layer certainly predominates. 
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3. Criteria for the Rupture of a Liquid Film on a Boiler Tube Surface 

It can be shown analytically and experimentally that there are significant 

differences between surface vaporization or wetted wall boiling, film boiling droplet 

14 

or rivulet flow and linear slug flow boiling . A major question that arises is, which 
of these modes of boiling heat transfer is in operation for a given set of parameters. 
One approach to this difficult area is to review a number of possible criteria that 
assist in defining conditions under which a liquid film that is in complete contact in 
a boiler tube (under surface vaporization conditions) ruptures. A number of obvious 
criteria are listed below and discussed briefly. 

Consider a perfectly smooth boiler tube surface upon which a liquid is 
vaporizing. It is postulated that the liquid is in direct contact with the surface 
and the boiling consists of surface vaporization or some volume boiling but no film 
boiling. If the contact angle between the liquid and the wall is zero degrees, then 
hypothetically, there will be a complete annular liquid film at the wall from qualities 
ranging from zero to 100 percent. If the contact angle (which is defined by the liquid- 
vapor, solid-vapor, solid-liquid surface forces) is greater than zero, the width of a 
liquid rivulet for that particular profile area defined by the contact angle and vapor 
quality must be compared to the inner circumference of the boiler tube. If the 
rivulet width is greater than the tube circumference then the liquid film would not 
rupture. If the rivulet width is less than the tube circumference, then the liquid 
film would be separated. 

Consider a rough boiler tube surface on which a liquid is vaporizing. 

Several possible processes can probably occur. For example, if the heights of 
the surface asperities are greater than the thickness of the liquid film, the contact 
area of the liquid film will be reduced. If the surface roughness is due to non- 
homogeneous structure such as oxide deposits or large metallic grains, changes 
in local contact angles can also cause local rupture of the liquid film, hence, reducing 
the contact area. Finally, as a result of surface pits or local areas of non-wetting, 
film boiling can occur at such points resulting in reduced superheat temperature 
differences at which surface vaporization ends and transition boiling begins. This 
change in boiling heat transfer mode is associated with a rupture of the liquid film. 
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4 . Critical Contact Angle at Film Rupture 

If a complete liquid annulus exists at a boiler tube wall and the liquid 
is in direct contact with it, the boiling conductance is high and increases 
with vapor quality. The heat transfer process involved is surface vaporization 
at the liquid-vapor interface or volume boiling . The important question of 
whether the liquid annulus is complete or not can hypothetically be answered 
by comparing the width of an axial rivulet relative to the boiler tube peri- 
meter. Further, the fraction of the boiler tube area covered with the 
direct-contact rivulet and the corresponding heat transfer are dependent 
upon the rivulet contact angle (defined by the liquid- vapor, solid- vapor 
and solid-liquid surface forces) and the vapor quality. 

Consider a liquid axial rivulet in a boiler tube. As the contact 

angle is decreased for a given liquid cross sectional area, the rivulet 

width increases and finally attains a value equal to the inner perimeter 

of the boiler tube. At this point the two edges of the rivulet meet and a 

complete liquid annulus is generated. The particular contact angle under 

this condition can be called the critical rupture angle (i.e. , if the angle 

is increased slightly, a hydrodynamic or surface disturbance could cause 

the annular film to rupture). The critical angle can be calculated in the 

following way. It is postulated that the rivulet profiles are controlled only 

by surface tension forces (all gravitational forces are presumed to be small) . 

Consequently, the rivulet cross sections are circular in shape. A second 

approximation involves the cross sectional area of a rivulet. It is postulated 

that the rivulet cross sectional area is the same as the value for a complete 

liquid annulus. From a previous Geoscience study, a liquid annulus boiling 

model has been developed that relates liquid film thickness at the boiler 

tube wall (or cross sectional area) and the heat transfer conductance to 
14 

vapor quality ; the results of that work are used in the present study. 

Consider the idealized rivulet system in Figure 96. The width of the 
rivulet at the critical contact angle is equal to the perimeter of the boiler tube, 
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(186) 


2r sin a 
s c 


2?rr t 


where 

r radius of curvature of rivulet surface 
s 

a , critical contact angle at film rupture 
c 

r^ , inside radius of the boiler tube. 

The rivulet cross sectional area is, 


ITT 


2 a 

c l 

— — - — (2r sin a 
27 r 2 v S c 


r cos a ) = 
s c 


2 7rr^. 6 


(187) 


where 6 is the thickness of the annular liquid film (a contact angle of zero or 

complete wetting). The solution of the two algebraic Equations (186) and (187) 

two unknowns (a and r) yields a relation for the critical contact angle, 
c 

26 2 1 . 

(sin a ) = (a — sin 2a ) (loo) 

iK t c v c 2 c' 


A calculation of the critical contact angle, based on this model and on the 
annular flow analysis of Reference 14, is presented at the end of the next 
section (page 213). 

5 . Rivulet Widths (Portion of Boiler Wall Covered by a Direct Contact 
Rivulet) 

Consider the case where the rivulet contact angle lies in the region 

0^<a< 90° (see Figure 97). If, as in the preceeding section, it is postulated 
that the liquid rivulet cross sectional area is equivalent to the cross sectional 
area for the case of complete wetting (annular liquid film) at a specific vapor 

q - aIity ’ 2 2a l 

7 ir - — (2r sina r cosa) = 2irr,6 (189) 

s 2 tt 2 ' s s t 


where r is the unknown and a is known for a given liquid, surface and tempera- 
s 

ture level. The fraction of the boiler tube area covered by a rivulet is, 

2r sina 


F 


s 

2 ^ 


(190) 
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Upon solving for r in Equation (189) and substituting it into Equation (190) there 
results 


F 


sing 
7 rr 




2 nr ^ 6 

a - — sin 2a 
2 


(191) 


A similar analysis for the case 90 < a < 180° yields 


F = sin (g ~ 


in (jr - a) . / 2rr t 6 
— 1 / — 
t T a+— si 


sin (jt - a) 


(192) 


As an illustrative example of the use of this analysis, consider a viscous- 

viscous annular flow pressure drop and heat transfer case from the analysis of 

Reference 14. Results for 6 and h ... (based on surface vaporization at 

an. iilm 

the film- vapor interface) obtained as a function of vapor quality , are 
listed in Table II. It should be noted that for the postulated viscous— viscous 
flow model, the working fluid flow rate does not enter the solution and is 
hence not needed. 


vapor quality, x 


0.1 


0.2 


0.3 


0.4 


Table II. 

5, ft h ... , Btu/hr ft 2 °F 

an. film 


0.0052 

0.0037 

0.0029 

0.0024 


3,270 

4,590 

5,860 

7,090 


Figure 98 shows the results of critical rupture contact angle ct^ as a 
function of vapor quality for the conditions of the illustrative example. For 
each vapor quality listed in Table II, a critical contact angle was calculated from 

Equation (188). If actual contact angles for the system fall below the critical values, 
the annulus film would not rupture; if the actual contact angles fall above the 
critical values, film rupture occurs. Note that critical contact angles decrease 
with increasing vapor quality. Thus for invariant actual contact angles, one 
would expect the annular film to rupture at higher vapor qualities. 

Figure 99 presents the function F versus vapor quality for contact 
angles greater than the critical angles for rupture. These values were cal- 
culated from Equations (191) and (192). 
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Figure 98. a c versus vapor quality for a potassium system 
(illustrative example). 
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Figure 99. Function F versus vapor quali 
(illustrative example). 
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6. Effect of Boiler Wall Thermal Conductivity on the End of the Wetted 

Film in a Boiler Tube 

Many experiments have shown that the onset of film boiling occurs when 
the superheat of the wall surface above the liquid saturation temperature reaches 
a certain temperature level. 

It was shown in Section D of this report that local heat fluxes and wall tem- 
peratures are much higher than the average values near discontinuities in the 
inner wall conductances. One might expect that when the surface temperature 
under the edge of the liquid film reaches a critical level local "vapor blanketing" 
will occur, and the end of the liquid film will move upstream. 

It might appear that in a uniform flux boiler the temperature under the 
edge would not change as the film edge moved upstream, and therefore the 
wall would dewet all the way to the beginning of the boiler. There is at least 
one plausible reason why this need not be the case. As stated in Section E.l, 
the unboiled liquid not carried in the liquid film is entrained in the vapor. The 
rate at which the entrained liquid is deposited on the wall is dependent on its 
concentration in the vapor, and the closer to the beginning of the boiler the 
film ends the higher the droplet concentration and the deposition rate must 
be. The impinging droplets do not wet the wall according to the model under 
discussion but they still contribute significantly to the wall cooling. The ratio 
of the 'wetted wall' heat transfer conductance to the 'dry wall' conductance 
would therefore be expected to decrease as the end of the film moves toward 
the boiler inlet. As is shown in Figure 69, the decrease in this ratio causes 
a lowering of the wall temperature at the edge of the film, so that at some 
point the edge becomes stable and ceases to move. Given a fixed value of 
wall to saturated fluid for which the onset of film boiling occurs, the maxi- 
mum heat flux for which the end of the annular liquid layer is stable is given 

by 
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(193) 


(q/A) 


max 


( T 1 - T s ) 
AT. 


max 


(lc/6) 


where (T -T ) is the temperature difference at which "vapor blanketing" occurs. 

1 s max 

Consider the application of this criterion to a potassium boiler. Suppose that for a 
specific case good estimates for , h^, and k/ <J have been made and that the values are 

h = 2 x 10 4 BTU/ft 2 hr *F 

h = 2 x 10 3 BTU/ft 2 hr °F 

k/<5 = 1 x 10 4 BTU/ft 2 hr °F 


Suppose also that the maximum superheat for which a wet wall may be maintained 

(T -T ) , is 50 °F. From the specified conditions and by use of Figure 69, 

x s max v 

AT X « 1.1 

so that 

(q/A) = 31 x 10 4 BTU/ft 2 hr 

max 

The specific case chosen corresponds to the temperature distribution given 
by Figure 67. 

The isotherm labeled 1 . 0 has a temperature (q/ A) / (k/ 6 ) above that of 
the inside wall at x = - °°. For the maximum heat flux this temperature differ- 
ence is 31 °F. This is also the temperature drop across the wall at a large 
distance from the discontinuity in either direction. 

Several conclusions may be made about the onset of film boiling as 
described by this model. One is that a thick walled or high conductivity boiler 
decreases the likelihood that the annular liquid layer will be removed by film 
boiling. Another conclusion is that film stability studies, such as that of Refer- 
ence 20 f i n which the attempt is made to describe analytically the conditions 
under which dry or wet spots will be maintained on a boiler wall, are not corn- 
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plete unless they include the effects of wall heat conduction. Still a third 
conclusion is that very rapid regime changes in a boiler tube might be 
interpreted as fairly gradual ones under some circumstances, if inside wall 
conductances are inferred from outside wall temperatures and fluxes. An 
example of such an interpretation may well occur in the boiling Freon data of 
Reference 21. Their data, obtained with a (high conductivity) copper tube, 
show a change from the maximum near the end of the annular film to the 
minimum in the fog flow region which is quite gradual, while the data for 
Freon given in Reference 22, obtained with an aid of flow visualization in a 
(low conductivity) glass test section show a discontinuous transition between 
the two regimes. A comparison of these data is shown in Figure 100. 

If instead of an annular liquid layer wetted liquid rivulets occur on the 
inside of the tube, as described in Section C. , lower outside heat flux would 
probably be required to bring on film boiling, since the variations in the local 
inside temperature differences are larger for rivulet flows than for discontinuous 
annular flows. 
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FREON 22 
FROM REF 21 
5.5 ATM. 


PARAMETER: G IN I0 4 lb /ft 2 hr 





F. Heat Transfer in Boiler Tubes 


1. Direct Contact Rivulet Heat Transfer 

The mean boiling heat transfer conductance in a boiler tube with a 
direct -contact axial rivulet can be expressed as 


h b - "l F + h v d - 


( 194 ) 


where F is the fraction of the tube wall covered by liquid, h^ is the liquid rivulet 
vaporization conductance, and h^ is the convective conductance for the dry areas 
in the boiler tube. The rivulet conductance can be defined in terms of the mean 
thickness of the rivulet, i and the rivulet thermal conductivity, . 


h 


(195) 


riv 


An alternate way in which to express the rivulet conductance is in teriAs of the 
conductance for a complete annular film and the function F, which is given on 
page 213. Thus, k 




l 


U riv 

2 r sin a 
s 

2 7rr. 


= h F 

an. film 


Combination of Equations (196) and (194) yields 


h = h ... F + h (1 - F) 
b an. film v 


(196) 


(197) 


The elementary heat transfer analysis presented above for a smooth 
boiler tube may be applied to an illustrative example. Consider the case of 
potassium forced flow boiling in a 3/4 inch ID tube, under the same conditions 
as were used to calculate the conditions of Table II on page 213. 
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Figure 101 shows how the boiling conductance function is degraded as 

the liquid annular film is ruptured and the rivulet contact area F is decreased. 

There is evidence in the case of a number of mercury experiments that this 

23 

phenomenon plays a major role in controlling the boiling heat transfer 
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VAPOR QUALITY, X 


Figure 101. Boiling conductances versus X and a for a potassium system 
(illustrative example) . 






2 . Effect of Surface Tension and Gravity on Rivulet Heat Transfer 

In the preceeding section, wherein surface vaporization heat transfer 
from liquid rivulets making direct contact with the boiler tube wall was considered, 
it was postulated that only surface tension forces were in operation to control 
the rivulet profiles. In the case of a high surface tension liquid such as mercury, 
this postulate seems to be reasonable. Under high velocity rotational flow 
conditions, however, a centrifugal force field is acting on a rivulet or droplet 
moving helically along the boiler tube surface. Therefore, it is desirable to 
calculate under what conditions the centrifugal acceleration field will significantly 
change the rivulet profile from what it would be if only surface tension were 
effective. The effect of gravitational forces (or centrifugal forces due to rotation) 
previously excluded are added in the rivulet profile analysis given below. 

From classical physics, the differential equation that describes the force 
balance between surface tension and gravitational forces for a rivulet (two dimen- 
sional system shown in Figure 102) is 




2-i3/2 
+ P )\ 



where 


P 


di 

dx 
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Figure 102. Rivulet coordinate system. 
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where 


y, vertical rivulet surface coordinate (from plane of zero pressure) 
x, horizontal rivulet surface coordinate (from rivulet center) 

c 

a , surface tension 
Pg , density 

g, acceleration of gravity 

Two boundary conditions define the system, namely, 

0 at the rivulet center (x = 0, y = y.) 

cot a at the rivulet edge (x =» a, y = y ) 

cl 

a , complement of the contact angle 

a, half width of rivulet 

y , vertical coordinate at edge of rivulet 
a 

The first integration of Equation (198) together with the first boundary condition 
yields 


dy 

dx 


and 


dx 



dx 


2 2 
c - 2 (y 


7 


72 2; 27 

(y - y 0 ) - (y - y D ) 


* 


dy 


(199) 


The final integration is facilitated by making the change of variable 


2 

y - 



2 2 
c cos $ 


( 200 ) 


Upon substituting Equation (200) into Equation (199) together with the definition, 

2 2 2 2 
k =0 /(y + c ), there results, 

o 


dx 


2 2 .2 
c - 2c sin <f> 

Y Vl - sin Z 4, 


d(J> 
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( 201 ) 


Because cp = 7r/2 when y = y or x = 0, the limits of integration are x to 0 and cp to n/2. 

o 

The resulting solution can be shown to be 


x = 


(i - y)( K(k) " F <M)) " ( E(k) “ E 


( 202 ) 


where the functions K(k), F(k,<£), E(k) and E(k,0) are the classical elliptic integrals. 

It is noted that both k and <p in Equation (202) are dependent upon. It is necessary to 

use the second boundary condition to complete the solution. Specifically, the value 

of <p at x = a, y = y , is given by 
a 


cos <p 


1 



(203) 


But since p = cot a when y = y , it can be shown that 

a 

2 

Vy - y 

a o_ _ 1 - sin a (204) 

c 2 

Thus, Equation (203) can be expressed as 


cos <p. 




- sin a 


Also Equation (202) becomes 


(205) 


a 


(k ' f )[ K<k) ' 



E (k,0 x ) 


(206) 


The solution for the profile of a rivulet is then obtained as follows. First, the value 
of k is determined from Equation (206) , after substitution of the constants <f> (from 
Equation (205) and c. The value for y^ is next determined from the definition of 
the term k. Finally, the values for x and y are given in terms of parameter (p by 
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Equation (202) and (200) respectively. The vertical coordinate at the edge of the 
rivulet, y , is given by Equation (204). 

d 

The cross sectional profile area of the rivulet is obtained for the complete 
solution outlined above by the relationship 

A = 2a y - 
a 

In order to compare this result to the analysis for which gravity was 
neglected, the width of a rivulet with a circular cross section (which neglects 
gravitational effects) has to be determined. That result is then compared to the 
exact value, 2a. 

The solution presented in this section has been evaluated for a number 
of specific mercury flow cases. Figure 103 shows a plot of the profile of a 
1.27 cm wide mercury rivulet under a one earth gravity field condition and a 
contact angle of ten degrees. The area under the curve was integrated and the 
result used to obtain a rivulet profile with a circular shape (no gravitational forces 
considered) with a contact angle of ten degrees. A second evaluation (Figure 104) 
pertains to a condition where a one hundred fold earth gravity field was in 
operation on a 1.27 cm wide mercury rivulet with a ten degree contact angle. 

Note the drastic changes in the rivulet profiles for these specific conditions. 
These results, which are particularly applicable to rivulet flow on the inside of a 
boiler tube wall with linear or helical flow, show that as one might expect, 
significant flattening of the rivulets occurs under high gravitational fields. At 
one earth gravity field, the ratio of rivulet width to the corresponding value for a 
zero gravitational field is 1.17. At one hundred earth gravity fields, the ratio of 
rivulet width to the corresponding value for a zero gravity field is 1.78. These 
results can be used to refine the rivulet heat transfer model presented in Section 
E.5. 
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3. Heat Transfer with Film Boiling Entrained Droplets in Helical Flow 

A recapitulation of Geoscience ' s boiling potassium and mercury mea- 
surements in a boiler tube with a twisted divider indicated that the boiling con- 
ductances generally decreased from predictable entrance values (at a vapor 
quality of zero percent) to predictable helical vapor flow exit values (at a vapor 
quality of one hundred percent); these changes in conductance were found to 
be two to three orders of magnitude (Reference 14). Some typical data are shown 
in Figure 105. 

The boiling heat transfer in the vapor quality region, 0. 5 <x < 1. 0, may 
be examined on the basis of two models previously proposed by Geoscience. 

One model postulates the presence of droplets or rivulets at the duct wall with 
a helically flowing gaseous core. A second model presumes that the drops are 
very small and distributed in the form of fog within the helically flowing super- 
heated gaseous core. Both models require that the liquid phase distribution 
and droplet size be known. On the basis of some earlier experimental measure- 
ments with mercury and air, some data on the droplet number per diameter 
length of boiler tube are available for helical flow; information on fog particle 
size and distribution is only qualitative, however. 

Consider the case where droplets and vapor are flowing helically in a 
boiler tube. Heat transfer from the wall to the working fluid consists of both 
gaseous forced convection and droplet vaporization. Because the mechanisms 
are operating in parallel, the corresponding conductances can be added, 

k = h +. h (208) 

d conv vap ' 

where h^, total boiling conductance 

h , convective conductance 
conv 

, droplet vaporization conductance 
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Figure 105 . Some typical potassium and mercury boiling conductance 
data in the high vapor quality region (helical flow). 
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that 


From the forced flow drop vaporization model derived earlier * it was shown 
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Upon substituting Equations (210), (211), and (212) into (209) and simplifying, 
there results, 
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vap 
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6/5 

X 


The quantity n varies with tube vapor quality because the vapor volume increases 
thereby increasing the distance between the droplets. It can be shown that, 

n x 

o o 

n = -7— (214) 
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where n , number of droplets per diameter of tube length at vapor quality, x 
o o 

x , vapor quality near tube entrance 
o 

Upon substitution of Equation (214) into (213) and simplifying Equation (212) for 

x > 0.01 and w, /w ~ 4000, one obtains, 

£ v 


h 

vap 


( 215 ) 
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The convective conductance term in Equation (208) has previously been 
expressed as 

^conv - ^forced + ^free ( 216 ) 

The conductances can be expressed in terms of the Nusselt modulus (Reference 14) 
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(218) 


Substitution of Equations (215), (216), (217) and (218) into Equation (208) gives an equation for 
the total boiling heat transfer conductance. Equation (208) has been compared in Figure 106 
to some experimental boiling potassium data obtained previously at Geoscience; the 
defining conditions for the forced, helical flow boiling are noted in the figure. The 


235 



Oi-er 



PARAMETERS: 
n =10; X 0 = O.I 
L = 850 BTU/LB 
/° d = 9.97 x I0‘ 8 LB-HR 2 /FT 4 

P - 5.58x10“" LB-HR Z /FT 4 

k v = 0.009 BTU/HR.FT °F 
W t = 13.3 LBS/HR 
At = 500°F 
4.23 


EXPERIMENTAL CHECK POINTS 


SYMBOL 

LBS /HR 

■ 

15.4 

• 

17. 1 

▲ 

8.8 


D = 1.58 xlO‘ z FT 
W v =0.0232 LBS/FT 3 

Wy 

— = 5.6 x lO' 4 





predicted and experimental data have similar profiles and magnitudes. The 

calculations were made for the case, n = 10 and x =0.1; these values were 

o o 

estimated on the basis of adiabatic mercury-air flow experiments in a glass 

, i 24 

tube 

4 1 Effects of Droplet Bouncing on Film Boiling in Helical Flow 

Previous Geoscience studies of the film boiling droplet vaporization 

regime indicated that a tiny droplet that was vaporizing in the steady state regime 

could, at times, suddenly begin to oscillate (to a height equal to, or greater than, 

24 

ten times the diameter of the droplet) • When this process occurs (for reasons 
not completely clear) , the corresponding vaporization lifetimes are found to be 
as much as fifty percent greater than those predicted for steady state film boil- 
ing. It is noted that the steady state predictions were usually in good agreement 
with steady state experimental droplet film boiling. Since the analysis in the 
preceeding section is based upon an analysis of steady state droplet vaporization, 
it is important to assess the effect of bouncing upon the helical flow film boiling 
analysis. 

The bouncing droplet dynamics in the film boiling state have previously 

24 25 

been discussed in Geoscience reports ’ . In the present section, two simple 

analyses have been made of the effect of a rotational force field on a droplet that 
leaves the surface of a boiler tube with a given, invariant velocity. The solutions 
obtained can be used to calculate the reduction in the height to which a bouncing 
droplet will rise for constant escape velocity at the wall. The results of the two 
analyses are outlined below, together with several restricting postulates. 

This complicated dynamics problem was idealized to the case of a simple 
bouncing sphere having a coefficient of restitution equal to unity and an invariant 
initial velocity (not related to the gravitational force field strength). It is recognized 
that these conditions may be too limiting for the purpose of drawing general conclusions 
at this time. 

First the classical situation of a uniform gravitational force field (with elevation) 
is considered. The highest point in the trajectory, y^, for this condition is 



y i " 2g (219) 
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or 






1 

t/s e 


where 

y , invariant initial velocity at y = 0 
o 

g , acceleration of the earth's gravity field 
€ 

g, acceleration of gravity for system 


( 220 ) 


A graph of Equation (220) is shown in Figure 107; note that for the postulated conditions, 
the absolute and dimensionless highest point in the trajectory varies inversely with 
the gravity field. Therefore, in a boiler tube with high velocity rotational flow, the 
high-amplitude bouncing droplet mode observed in the one-gravity droplet vaporiza- 
tion apparatus is not expected to be as prevalent. 

It is likely that the gravitational field independent of height is too restrictive 
for the present purposes, so a second analysis was made for the case of a variable 
gravitational force field as found in a boiler tube having helical flow . The general 
force equation on the droplet is 


w •• 

— tt y = - W 
gy 


( 221 ) 


where 


W, weight of the droplet 

g(y) , variable acceleration field 

y, second derivative of y with respect to time 

t, time 


Upon expressing the first derivative of y with respect to time as the variable p, 
Equation (221) is transformed to 

P ^ = - g(y) (222) 

This equation may be integrated twice to yield 
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(223) 


I.T. 


dy 


y]~ 2 f g(y) dy 


+ 2C 


/ d 


dt + C 


2 

From the condition y = at t = 0, it can be shown that the constant = l/2 (y Q ) . 
From the condition y = 0 at t = 0, it can be shown that the constant = 0. From 
the first derivative of y with respect to t, the amplitude of the trajectory, y , can be 
found from the equation 


*1 



g(y) dy = y 


(224) 


For helical flow in a boiler tube the tangential velocity, v^ can be approximated by 


v 

t 


= \ <r t 


- y) 


(225) 


where 

b^, a constant dependent on the axial flow rate 

r t , boiler tube radius 

y, radial distance from the tube wall 


The centrifugal acceleration term is considered to control the acceleration g(y) for 
most practical cases; thus, 


g(y) = 




\ 2 ( r t - y) 


(226) 


Substitution of Equation (226) into Equation (224) leads to 



r t y i 



( 227 ) 


The solution of this quadratic equation for the region 0 < y^/ r < 1 is , 
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t 


1 



(228) 


where the quantity, y / v t » is the ratio of vertical to tangential velocities at the 

° o 

boiler tube wall. A graph of Equation (228) is shown in Figure 108. Note that 
the trajectory amplitude of a droplet decreases as the tangential velocity 
(acceleration field) increases . This is in agreement with the result of the 
analysis for a uniform force field, the results for which are given in the 
preceeding figure. 
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5 . Heat Transfer with Deposition and Vaporization of Entrained Droplets 

One possible mechanism for the heat transfer behavior in the "transition 
boiling" region which occurs downstream of the end of the annular film, is 
spray cooling of the wall by wetting evaporation of droplets deposited from the 
entrained flow. 

An analysis has been presented for the rate of deposition of entrained 
droplets in Section B.6. If one postulates that the wall is locally 
wetted by the impinging droplets, and that they are completely vaporized on the 
wall, then a computation may be made of the effect the droplet deposition will 
have on the mean temperature of the vapor flow in the tube. The temperature of 
interest here is the mixed mean temperature of the vapor phase only. The droplets 
entrained in the vapor flow are at saturation. The idealizations upon which the 
analysis is based are: 

(1) The annular film has ended 

(2) There are no rivulets; all of the remaining liquid is entrained in 
the vapor flow 

(3) Deposition follows the relationships given in Section B.6 

(4) The droplets wet the wall and vaporize completely 

(5) The wall is subjected to a uniform heat flux 

(6) The vapor temperature is initially at saturation 

(7) No significant vaporization occurs in the entrained flow. 


A simple energy balance based on these idealizations leads to the following 


equation for the mixed - mean vapor temperature: 



(229) 
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where 


T = 


T " 


L 


c 

P 


T 

s 


q/A n D£ 

k = ~ 

2 L w 

eo 

L = latent heat of vaporization, c = vapor specific heat, and s and w are the 

p o eo 

normalized boiler length and the entrained liquid flow at the end of the annular film. 
The other symbols are defined in section B.6. 


Inspection of Equation (229) shows that in the absence of deposition, 
the mixed mean vapor temperature rises linearly, consistent with the ideali- 
zations. Physically the smallest value of T which may occur is zero, and this 
value will occur if the deposition (and vaporization) rate is large enough to re- 
quire all of the heat provided to the wall. Mathematically, the value of T can 
be negative, since the mathematical relationships can specify a deposition rate 
larger than that which can be completely vaporized. In this latter case, a li- 
quid film would form on the wall, which is a violation of the postulates. 

Two cases have been computed, one for Freon 12 and one for potas- 
sium. These cases correspond to the ones chosen for the deposition rate cal- 
culations presented in section B.6. The results are presented in Figures 109 
and 110. 

The negative temperature in Figure 109 shows that according to the 
model a liquid film re-forms on the wall. Since the case chosen corresponds 
to a "burnout" point presented by Kirby (Reference 12) and therefore accord- 
ing to experiment the annular film has in fact ended, it is clear that the depo- 
sition rate is too high, (by at least a factor of four). . The most likely reason 
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Figure 110. Vapor temperature rise downstream of end of annular flow. 
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for the overestimation of deposition rate is that the diffusion conductance given 
by Equation (80) of this report is too large for the present case. Equation (80) 
is appropriate for very small particles, which follow the turbulent eddies of 
the flow in the tube. The droplets known to exist in a boiling Freon flow are 
too large to follow the local eddy motion and do not, therefore, diffuse as rapid- 
ly. In addition, the flow in an actual boiler is accelerating, and has mass 
addition next to the wall, so that the actual turbulence level is not that of the fully 
established pipe flow which Equation (80) represents. It may be concluded that 
although a large effect of spray cooling on transition boiling heat transfer may 
possbily be expected in a Freon boiler, it is overestimated by the present 
analysis* 

Figure 110, presenting the calculation results for potassium, shows a negative 

normalized temperature also, indicating that a liquid film is forming faster than 

it is being evaporated away. It may be seen in the figure that the mixed-mean 

vapor temperature begins to rise at s - s =0.7, which is approximately 92 

o 

diameters downstream of the point s^. At the point where the mixed-mean vapor 
temperature begins to rise the vapor quality is about 81%. 

Although the analysis probably overestimates the droplet deposition rates 
for both Freon and potassium, it appears likely that, in some circumstances, 
deposition rates could be high enough to significantly affect the conductance down- 
stream from the end of the annular film. Droplet deposition probably plays a role 
in some cases of transition boiling. 

The postulate was made in the analysis that the liquid wets the boiler 
wall. It is well known that wetting will occur only if the wall temperatures are 
not too high above the saturation temperature. Under most conditions of uniform 
wall heat flux, the wall temperature will rise rather rapidly after the end of the 
annular film and somewhere downstream the droplets will no longer wet the wall, 
but will instead film boil. 

At this point the present analysis will no longer apply. The calculation 
of boiler wall temperature rise downstream of the end of the annular film may 
also be performed by use of the approximation that the vapor forced convection and 
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Reference 14, the rivulet flow analyses of section F.2, the criterion for onset 
of nucleate boiling (section E.2), suitably modified for appropriate nucleating 
conditions, and the varying wall temperature information of section D. 

The use of simple models and their combinations to make performance 
predictions appears to be the most promising method of avoiding a cut-and-try 
experimental program for each new requirement. Test programs are essential, 
but with an initial performance prediction, and careful analysis of the test results, 
using analysis as a guide, the test program may hopefully be short, and rational. 
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IV. NOMENCLATURE 


A heat transfer area 

A liquid droplet cross-sectional area 

d 

A tube wall inside area 

c 

A, tube wall inside area wetted by rivulet 

i 

A tube wall outside area 

o 

A rivulet cross-sectional area 

r 

A_ tube cross-sectional area 

T 

a 1) a constant; 2) rivulet half width 

a radial acceleration 

r 

B 1) a constant; 2) a parameter 

b a constant 

b t a constant 

C, C C 2 , Cg. C 4 constants 
drag coefficients 

C D wire drag coefficient modified for boundary layer effects 
m 

drag coefficient for a circular cylinder adjacent to a solid wall 
w 

0 ^ vapor specific heat (constant pressure) 

C g vapor stream droplet concentration 
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C droplet concentration near tube wall 
w 

D tube inside diameter 

D tube equivalent diameter for linear flow 

e 

D tube equivalent diameter for helical flow 
h 

D wire coil inside diameter 

m 

^ droplet (or particle) diameter 

tube equivalent diameter for tangential flow (normal to tube axis) 

2) molecular diffusion coefficient 

d helical wire diameter 

w 


E (K) 

E (k,0) 


elliptic integrals 


F 1) force; 2) fraction of boiler tube inside area covered by rivulets 

F (k, <f> ) elliptic integral 

F drag force 

F wall friction force measured normal to tube axis 

F gravitational force 

g " 

F centrifugal force due to rotation of droplet around tube axis 
H 

F force exerted on fluid by helical wire measured normal to helical wire 
n 

F 1) forces restraining secondary flows; 2) radial shear forces on film due to 
secondary flows 

F force driving secondary flows 

s 
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wt 


wire drag force measured normal to tube axis 
F^ vapor drag force in x direction (droplets in helical flow) 

F vapor drag force in y direction (droplets in helical flow) 

y 

G mass velocity (of both phases) in boiler tube 

g local acceleration field 

g (y) variable local acceleration field 

g proportionality constant in Newton's second law 

c 

g Earth' s gravitational field 


H normalized conductance, defined by H = h/h 


x = 100% 


heat transfer conductance 


h apparent conductance on inside tube wall, observed from outside 
a 

h heat transfer conductance with a complete annular film 

an. film r 

h^ boiling heat transfer conductance 

h convective conductance 

conv 

h , a mass diffusion conductance 
d 

h. , single phase forced convection conductance 
forced 

h^ee convective conductance due to buoyancy (natural convection) 

h^ convective conductance on hot fluid side of heat exchanger surface 

h„ wetted wall conductance 

i 

h^ dry wall vapor conductance 

hyap conductance due to droplet vaporization 
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I 1) rotational moment of inertia; 2) an integral function 

K (k) an elliptic integral 
k, k^, k^ parameters, constants 
k^ liquid thermal conductivity 

k^ vapor thermal conductivity 

k wall thermal conductivity 

L latent heat of vaporization 

a characteristic length 

A length 

A axial tube length 

a 

A, helical path length of vapor flow 
h 

A maximum altitude of droplet ejected from a horizontal liquid surface 

max 

A axial station at which droplet motion begins 

A characteristic lengths 

i transverse distance of droplet travel 

A length of helical wire 

w 

M applied torque 

m droplet mass 

Nf, N , N , parameters 
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Nusselt modulus 


N 

u 

N Nusselt modulus for a linear turbulent single phase vapor flow 

U lin.turb. vap 

N Nusselt modulus of forced convection 

U forced 

N Nusselt modulus of free convection 

U free 

n number of droplets per unit length of boiler tube , number of vanes on helical divider 

n number of droplets in an axial length equal to one tube diameter 

d 

n number of droplets per unit length at quality X 
o o 

n^, n^ parameters 


P 

r 

P 

P c 

P f 

P s 

P w 

P wm 

q/A 

(q/A) 

(Q/A) 

R. 


Prandtl modulus 

1) helical wire or vane pitch; 2) a point in space; 

pressure at centerbody surface 

pressure drop due to skin friction 

saturation pressure 

pressure at tube inside wall 

maximum pressure at tube inside surface 

heat flux 

"critical" or "burnout" heat flux 
B 

uniform outside heat flux 
o 

linear flow Reynolds modulus 
inside tube wall radius 
outside tube wall radius 


3) 


dx 


4) 


dy 

dt 


pressure 
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r radius measured from tube centerline 

r , r radius of inner channel wall or centerbody 

o 1 

r radius of inside surface of outer channel wall 
2 

r^ radius of centerbody 

r^ radial location of liquid-vapor interface 

r radius of curvature of rivulet surface 
s 

r radius of boiler tube (inside) 

Sc Schmidt modulus 

s normalized axial length 

normalized axial length of which annular film ends 

T absolute temperature 

a 

T mixed mean temperature of hot fluid 
h 

T outside wall temperature 

o 

T g saturation temperature 

mixed mean vapor temperature 
wall temperature 

T normalized mixed mean vapor temperature 

t time 

t time at which droplet reaches maximum altitude 

£max 

t normalized time for axial motion 

a 

t normalized time for radial motion 

r 

\l 80 ° normalized time at which droplet has traveled half way around tube 
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droplet volume 

V normalized particle velocity 

P 

V g normalized slip velocity 

normalized vapor velocity 
v velocity 

v q initial velocity of droplet ejected from liquid layer 

v^ mean axial velocity 

v critical velocity for shear entrainment 

e 

v, helical velocity at tube surface 

h 

V max maximum axial fl° w velocity in a tube with a velocity profile 

v^ velocity component normal to helical wire 

v droplet (or particle) velocity 

P 

v 1) a characteristic velocity; 2) velocity of liquid film on vane, measured in 
radial direction 

v 1) slip velocity; 2) secondary flow velocity 

s 

v^ velocity component normal to tube axis 

v initial transverse droplet velocity 

v^ velocity component normal to tube axis in fully established flow 

v transverse droplet velocity 

v vapor velocity 

v x component of vapor velocity 

v y component of vapor velocity 

v velocity of liquid film on vane, measured in direction of constant radius 

X 
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W weight of droplet 

W t total mass flow rate 

w angular velocity 

w mass flow deposited onto tube wall 
d 

w entrained liquid mass flow rate 
e 

w initial entrained liquid mass flow. rate 

eo 

Wg weight density of liquid 

w^ net flow of liquid film in radial direction 

1) vapor phase flow rate; 2) weight density of vapor 


X 


X_ 


B 


X 


x 


X 

o 


X 

r 


x 

e 


vapor quality 
quality at "burnout" 

quality at which number of droplets is known 
a space coordinate 

coordinate of droplet motion initiation 
half width of rivulet 

distance from rivulet center to wall station where temperature gradient is zero 


Y„ void fraction 

£ 

y a space coordinate 

y first derivative of y with respect to time 

y second derivative of y with respect to time 

y^ initial velocity of bouncing droplet at y = o 
y^ highest point in the trajectory of a bouncing droplet 
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z a space coordinate lying along the boiler tube axis 
z a normalized axial coordinate 

axial position of initial droplet motion 

c t 1) a parameter; 2) contact angle 

a , ot „ parameters 

.L Z 

ot^ critical contact angle at film rupture 

volumetric thermal expansion coefficient 
a a normalized gravitational force parameter 

Z a difference quantity 

Z p^ w pressure drop in a tube containing a helical wire 

Z T wall-to-saturation temperature difference 

ZT^ temperature drop across boiler tube wall 

6 1) boiler tube wall thickness; 2) liquid film thickness 

6 . mean thickness of liquid rivulet 

rxv 

€ eddy diffusivity 

£ friction factor 

4^ friction factor based on tube diameter 

friction factor based on wire coil inside diameter 
m 

6 angle measured about the tube axis 

k a parameter 
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X a reciprocal characteristic time 

p a parameter 

p gas phase viscosity 

p^ liquid viscosity 

p vapor viscosity 

p fluid density 

p gas phase density 

5 

P ^ liquid density 

P y vapor density 

<7 surface tension 

t a characteristic time 

t a characteristic time for axial droplet motion 

a 

7* r a characteristic time for radial droplet motion 

t shear stress 

s 

0 a function 

6 angular position in the curved channel 
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